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ABSTRACT 

We aim at characterising dense cores in the clustered environments associated with 
intermediate/high-mass star- forming regions. For this, we present an uniform anal- 
ysis of Very Large Array NH3 (1,1) and (2,2) observations towards a sample of 15 
intermediate/high-mass star- forming regions, where we identify a total of 73 cores, 
classify them as protostellar, quiescent starless, or perturbed starless, and derive some 
physical properties. The average sizes and ammonia column densities of the total 
sample are ~ 0.06 pc and ~ 10 15 cm -2 , respectively, with no significant differences 
between the starless and protostellar cores, while the linewidth and rotational temper- 
ature of quiescent starless cores are smaller, ~ 1.0 km s _1 and 16 K, than linewidths 
and temperatures of protostellar (~ 1.8 km s _1 and 21 K), and perturbed starless 
(~ 1.4 km s _1 and 19 K) cores. Such linewidths and temperatures for these quiescent 
starless cores in the surroundings of intermediate/high-mass stars are still significantly 
larger than the typical linewidths and rotational temperatures measured in starless 
cores of low-mass star- forming regions, implying an important non- thermal compo- 
nent. We confirm at high angular resolutions (spatial scales ~ 0.05 pc) the correla- 
tions previously found with single-dish telescopes (spatial scales > 0.1 pc) between 
the linewidth and the rotational temperature of the cores, as well as between the ro- 
tational temperature and the linewidth with respect to the bolometric luminosity. In 
addition, we find a correlation between the temperature of each core and the incident 
flux from the most massive star in the cluster, suggesting that the large temperatures 
measured in the starless cores of our sample could be due to heating from the nearby 
massive star. A simple virial equilibrium analysis seems to suggest a scenario of a 
self-similar, self-graviting, turbulent, virialised hierarchy of structures from clumps 
(~ 0.1-10 pc) to cores (~ 0.05 pc). A closer inspection of the dynamical state taking 
into account external pressure effects, reveal that relatively strong magnetic field sup- 
port may be needed to stabilise the cores, or that they are unstable and thus on the 
verge of collapse. 

Key words: ISM: molecules - radio lines: ISM - stars: formation - stars: early-type. 



1 INTRODUCTION 
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The initial conditions of the star formation process in clus- 
ters are still poorly understood. Studies have unveiled the 
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Table 1. Massive star-forming regions observed in NH3 with the Very Large Array, and few observational and derived parameters 



VLA 



Phase centre b 



Synthesised beam c 



Clumpfmd 6 



Region a 


(kpc) 


(£©) 


project 


R.A. 


Dec. 


HPBW 


P.A. 


RMS d 


P 


B 


S 


1. 00117+6412 


1.8 


1300 


AB1217 


00:14:27.725 


+64:28:46.17 


3.9 


x 3.5 


+46 


2.2 


66 


35 


10 


2. AFGL5142 


1.8 


2200 


AZ120 f 


05:30:48.024 


+33:47:54.44 


3.6 


x 2.6 


-14 


3.5 


203 


20 


10 


3. 05345+3157NE 


1.8 


630 


AF471 


05:37:52.400 


+32:00:06.00 


2.3 


x 2.3 


+92 


1.4 


39 


25 


5 


4. 05358+3543NE 


1.8 


3100 


AC733 


05:39:12.600 


+35:45:52.30 


5.8 


x 4.1 


+53 


5.3 


269 


40 


10 


5. 05373+2349 


1.6 


1100 


AF484 


05:40:24.500 


+23:50:55.00 


3.1 


x 2.2 


-73 


4.0 


116 


30 


10 


6. 19035+0641 


2.2 


7700 


AF386 


19:06:01.610 


+06:46:35.80 


4.0 


x 3.7 


-13 


2.6 


378 


28 


10 


7. 20081+3122 


2.5 


28200 


AF386 


20:10:09.050 


+31:31:35.20 


5.0 


x 4.4 


-56 


3.3 


602 


25 


5 


8. 20126+4104 


1.64 


8900 


AZ113 


20:14:26.053 


+41:13:31.49 


3.6 


x 3.2 


+10 


3.8 


873 


20 


10 


9. G75.78+0.34 


3.8 


96000 


AF386 


20:21:43.970 


+37:26:38.10 


3.7 


x 3.4 


-57 


3.8 


181 


30 


5 


10. 20293+3952 


2.0 


8000 


AZ120 


20:31:10.698 


+40:03:10.75 


6.9 


x 3.1 


+72 


3.0 


240 


40 


10 


11. 20343+4129 


1.4 


1500 


AS708 


20:36:07.301 


+41:39:57.20 


4.2 


x 3.2 


+11 


3.5 


224 


27 


10 


12. 22134+5834 


2.6 


11800 


AK558 


22:15:08.099 


+58:49:10.00 


3.8 


x 3.1 


+88 


1.5 


31 


40 


10 


13. 22172+5549N 


2.4 


830 


AF484 


22:19:08.600 


+56:05:02.00 


3.2 


x 3.0 


+17 


1.2 


61 


30 


10 


14. 22198+6336 


0.76 


340 


AS926 f 


22:21:26.764 


+63:51:37.89 


3.7 


x 3.1 


-41 


3.8 


61 


40 


10 


15. CepA 


0.7 


25000 


AF386 


22:56:17.870 


+62:01:48.60 


5.4 


x 4.9 


-29 


3.6 


332 


30 


6 



a Name of regions starting with numbers (e. g., 00117+6412) refer to the IRAS name (Neugeb auer et al.l 1984). 

b Phase centre coordinates (in Equatorial J2000.0) of the ammonia maps. Right ascension in units h:m:s and Declination in units d:m:s. 

c Synthesised beam (HPBW) in arcsec, and position angle (P.A.) in degrees. 

d RMS noise level, in units of mJy beam -1 , per channel (of 0.3 km s _1 or 0.6 km s -1 ; see Table [2}. 

e P: intensity peak (in mJy km s — 1 ) of the zero-order moment map. B (bottom) and S (step) parameters, in percentage of the peak, 

used in the Clumpfmd algorithm to extract the position of the cores (see Sect. 13. lj) . 

f AFGL5142 and 22198+6336 includes also data of the project AZ114. 



physical and chemical properties of starless isolated low- 
mass cores on the verge of gravitational collapse in nearby 
low-mass star- forming regions, showing that they have dense 
(n ~ 10 5 -10^ cm" 3 ) and cold (T ~ 10 K ) nuclei (e.g., 
iTafalla et aDl200l 12004 ISchnee et aD 1201(f ). and that the 
internal motions are thermally dominated, as demonstrated 
by their close-to-thermal line wi dths, even when ob served at 
low angular resolution (see lBergin fe Tafallal 12007, for a re- 
view). In the cold and dense nuclei of these cores, C-bearing 
mole cular species such as CO and CS are strongly depleted 
(e. g JCaselli etHll 20021 : pTaf alia et all2002h . while N-bearing 
species such as N2H + and NH3 maintain large abun dances 
in the gas phase (e. g. JCaselli etTaill 20021 : ICrapsi et alJl2005h . 

One major issue that has so far been poorly investi- 
gated is if, and how, the environment influences the physical 
and chemical properties of these pre-stellar cores. In clus- 
ters containing several forming stars in a small volume (of 
diametre ^0.1 pc), turbulence, relative motions, and inter- 
action with nearby forming ( proto-) stars can affect the les s 
evolved condensations (e. g., I Ward- T hompso n et al.ll2007l ). 
Such an interaction is expected to be much more important 
in clusters containing intermediate- /high-mass protostars or 
newly formed massive stars, given the typical energetic feed- 
back associated with the earliest stages of massive star for- 
mation (powerful outflows, strong winds, expanding HII re- 
gions) and the high pressure of the parental pc-scale clumps. 
There is also vigorous theoretical debate on how star forma- 
tion proceeds in clustered regions. Specifically: are the local 
kinematics of the gas dominated by feedback from protostel- 
lar outflows of already- forming, gene rally low- mass stars o r 
from feedback from high-mass stars (jNakamura &; Lill2007l )? 
Can a model of star formation starting from quiescent star- 
less cores, which has been successfully developed to explain 
observations of regions of isolated low-mass star formation, 



be applied to clustered regions (|McKee fe Tanll2003h ? To 
constrain theoretical models, it is crucial to measure the 
main physical and chemical properties of starless and star- 
forming cores in clustered environments. 

An observational effort in this direction has started, but 
it is mostly concentrated on nearby low-mass star-forming 
regions like Ophiuchus (e.g., Andre et al. 2007: Friesen et 
al. l20Q9h and Perseus (e.g.. lFoster et al.ll200^ due to the 
fact that the typical large distances (^ 1 kpc) of high-mass 
star-forming regions ma ke the study o f clust ered environ- 
ment svervchallenging. F oster et al] (120091) and Friesen 
et al. (2009) show that starless cores (studied at spatial 
scales ~ 0.01-0.05 pc) within low-mass star-forming clus- 
ters have typically higher kinetic temperatures (~ 13-14 K) 
than low-mass isolated cores (~ 10 K). On the other hand, 
the kinematics of these low-mass protoclusters seem to be 
dominated by thermal motions like in more isolated cores, 
even though the external environment is turbulent (Andre 
et al. [2007|). On the contrary, in the very few published 
studies of (proto-) clusters containing an intermediate- /high- 
mass forming star, the internal motions of starless cores are 
dominated by turb ulence (e. g.j Fontani et al. 2008, 2009 for 
IRA S 05345+3157: IWang et al. II 20081 for G28.34+0.06: Pillai 
et al.l201ll for G29.96-0.02 and G35.20-1.74: lFontani et al.l 
l2012bl for IRAS 20343+ 4129). 

By using ammonia, IPalau et al.l (|2Q07bh measured the 
temperature of the starless cores in the high-mass proto- 
cluster IRAS 20293+3952 at spatial scales of - 0.05 pc, 
through Very Large Array observations of the NH3 (1,1) 
and (2,2) inversion transitions. In fact, NH3 is the ideal 
thermometer for cold and dense gas because it is produced 
by the volatile molecular nitrogen, which is not expected 
to freeze out even in very cold and dense gas, and the 
NH3 (2,2) to (1,1) line ratio is sensitive to temperature (e. g., 
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Table 2. Configurations used, epochs, calibrators and spectral setup of the Very Large Array observational projects 



VLA 



Epoch 



Flux 



Gain calibrators 



Bandwidth Spectral resolution 



Project 


config. 


of observation 


calibrator 


(Bootstrapped fluxes, Jy) 


(MHz) 


(kHz) 


(km s- 1 ) 


AB1217 


D 


2007 


May 


3C286,3C48 


J0102±584 


3.85 ±0.06 


3.125 


48.8 


0.6 


AC733 


DnC 


2004 


Jun 


3C84 


J0530±135 


3.33 ±0.02 


3.125 


48.8 


0.6 


AF386 


D 


2001 


Jan 


3C286 


J2025±337 
J1849±005 
J2322±509 


2.50 ±0.04 
0.81 ±0.01 
0.72 ±0.01 


3.125 


48.8 


0.6 


AF471 


D 


2009 


Oct/Nov 


3C286,3C48 


J0555±398 


3.15 ±0.01 


3.125 


24.4 


0.3 


AF484 


D 


2009 


Oct/Nov 


3C48 


J0559±238 

J2148±611 


0.80 ±0.01 
0.63 ±0.01 


3.125 


24.4 


0.3 


AK558 


D 


2003 


May 


3C286 


J2148±611 


0.60 ±0.01 


3.125 


48.8 


0.6 


AS708 


C 


2001 


Jul 


3C286 


J2015±371 


2.34 ±0.04 


3.125 


48.8 


0.6 


AS926 


C 


2008 


Apr 


3C286,3C48 


J2146±611 


0.72 ±0.02 


3.125 


48.8 


0.6 


AZ113 


D 


1999 


May 


3C286,3C48 


J2013±370 


2.27 ±0.21 


3.125 


24.4 


0.3 


AZ114 


D 


1999 


Mar 


3C48 


J2230±697 
J0552±398 


0.43 ±0.01 
3.46 ±0.05 


3.125 


48.8 


0.6 


AZ120 


C 


2000 


Apr 


3C48 


J0552±398 


4.72 ±0.15 


3.125 


48.8 


0.6 




D 


2000 


Sep 


3C48 


J0552±398 
J2015±371 


3.38 ±0.04 
3.84 ±0.02 









IHo &; Towned ll983). Moreover, because other volatile species 
(like N2H + ) cannot be used to derive the tempera ture, am- 
monia is the unique tracer for this purpose. IPalau et al.l 
(2007b) found temperatures of 16 K for starless cores in 
IRAS 20293+3952, i.e., higher than those measured in iso- 
lated low-mass starless cores (or infrared dark clouds with 
no active star formation; e. g., Ragan et al.ll201ll ). but similar 
to the values measured in active clustered cores in Perseus 
and Ophiuchus. 

In this paper, we study a sample of 15 massive star- 
forming regions (see Table [J) observed with the Very Large 
Array (providing an average spatial resolution of ~ 0.05 pc) 
and analysed them in an uniform way. The main goal is to 
understand whether the results obtained in the few exam- 
ples shown above are also obtained when analysing a statis- 
tically significant sample. The regions were selected accord- 
ing to the following criteria: i) regions must be at distances 
< 3.5 kpc, to obtain spatial resolution < 0.05 pc, similar to 
the s cale of the cores (0.03-0.2 pc: e.g., iBergin &; Tafallal 
120071 ); ii) regions must have bolometric luminosities larger 
than rsj 300 Lq, which are the l uminosities wher e cluster- 
in g; seems to be important (e. g., iTesti et al.lll 999: Palau et 
al. l2013f ); and iii) regions must be still associated with im- 
portant amounts of gas and dust judging from single-dish 
observati ons of dust and molecular gas studies (e.g., Moli- 
nari et al. 1996; B euther et al.l 12002 ' 



2 OBSERVATIONS 

The Very Large Array (VLA[j) was used to observe the am- 
monia (J,K) — (1,1) and (2,2) inversion transitions toward 
15 intermediate/high-mass star-forming regions. These ob- 
servations are part of different observational projects carried 



1 The Very Large Array is operated by the National Radio As- 
tronomy Observatory (NRAO). The NRAO is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



out on different epochs (from 1999 to 2009), and with the ar- 
ray in compact configurations (C, DnC or D). In Table [1] we 
list the 15 massive star- forming regions, indicating the name 
of the VLA project and the coordinates of the phase centre 
for each region. More technical details of each observational 
project, such as calibrators and epochs, are provided in Ta- 
ble [2] The full- width at half-maximum of the primary beam 
at the frequency of the observations (~ 23.7 GHz) is approx- 
imately 2 arcmin. The spectral setup configuration used in 
each project (i.e., bandwidth and spectral resolution) are 
indicated in Table [2] The absolute flux scale was set by ob- 
serving the quasars 1331+305 (3C286), 0137+331 (3C48) or 
0319+415 (3C84), for which we adopted a flux of 2.41 Jy, 
1.05 Jy and 16 Jy, respectively. Gain-calibration was done 
against nearby quasars, which were typically observed at 
regular time intervals of ~ 5-10 minutes before and after 
each target observation. Flux calibration errors could be of 
up to 20-30%. The data reduction followed the VLA stan- 
dard guidelines for calibration of high-frequency data, using 
the NRAO package AIPS. Final images were produced with 
the task IMA GR of AIPS, with the robust parameter of 
iBriggsl ([1995) typically set to 5, corresponding to natural 
weighting. Two sources (AFGL5142 and 22198+6336) were 
observed in different projects with different angular resolu- 
tions. For these sources we combined the uv-d&ta, to improve 
the sensitivity and m;-coverage of the final images. The re- 
sulting synthesised beams and rms noise levels (per channel) 
are indicated in Table [T] 



3 RESULTS AND ANALYSIS 

3.1 Core identification and classification 

In Fig. [1] we present the zero-order moment (velocity- 
integrated intensity) maps of the NH3 (1,1) emission for the 
15 massive star-forming regions studied in this work, while 
in Fig. IB1I we present the first-order (intensity weighted 
mean vlsr) and second-order (intensity weighted mean ve- 
locity linewidth) moment maps. The zero-order moment 
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Figure 1. Colour image: NH3 (1,1) zero-order moment (integrated intensity) maps in units of mjy km s — 1 (horizontal colour-bar on 
the top of each panel). Coordinates in the x and y axis are right ascension and declination offsets in arcsec, with the (0,0) corresponding 
to the phase centre of the map (see Table ^). Black thick lines show the polygons for each core listed in Table [3] as obtained from the 
CLUMPFIND analysis (see Sect. [3]). The numbers inside the polygons identify the core (see Table EJ). The synthesised beams (listed in 
Table ^) are shown in the bottom-right corner of each panel. Star symbols in all panels correspond to centimetre continuum sources. 



maps were constructed by integrating the emission of all 
the hyperfine components, while for the first and second- 
order moment maps we considered only the main hyperfine 
component, except for 20126+4104 for which the inner satel- 
lite hyperfine component was used (in the zeroth, first and 
second order moments) to avoid the high opacities in the 
main line. The integrated intensity maps reveal ammonia 
emission in all the regions, coming from compact condensa- 
tions surrounded by faint and more extended emission. In 



the first-order moment maps (see Fig. lBl[h we identify large- 
scale velocity gradients (see e.g., AFGL5142, 20081+3122, 
G75. 78+0.34, 22172+5549N), not only associated with the 
compact condensations or cores but with the whole ammo- 
nia emission, probably revealing the large scale kinemat- 
ics of the cloud. To extract and identify the cores within 
each regi on we ran the two-d imensional CLUMPFIND al- 
gorithm d Williams et all 1 19941 ) on the NH 3 (1,1) zero-order 
moment map. This algorithm requires two input parameters: 
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Figure 1 - continued 



the 'bottom' level, corresponding to the lowest intensity to 
be included in the core, and the 'step', which determines 
the intensity increment req uired to differentiate o ne core 
from another core (see e.g., iMiettinen et al.ll2012i for fur- 
ther details). In our case, we used 'bottom' values of around 
20-40% of the peak intensity of the zero-order moment, and 
'steps' in the range 5-10% of the peak intensity (see last 
columns in Table [1]) . The output parameters defining each 
core returned by CLUMPFIND are the peak position, peak 
intensity and the flux density of the core. For each identi- 
fied core we defined a polygon at the level of the full width 
at half maximum (FWHM; see black thick lines in Fig. [1]). 
In Table [3] we list the offset position of each core (relative 
to the phase centre; listed in Table [1]) and the deconvolved 
FWHM size. A total of 73 cores were identified. 

The studied field of view (VLA primary beam at 1.3 cm) 
is ~ 2 arcmin, or 1.5 pc at a distance of 2.5 kpc. On average, 
we found of the order of 5-6 cores in each region, with di- 
ametres ~ 0.06 pc. We classify the identified cores as starless 
or protostellar (see Table 2} , based on their infrared proper- 
ties (from the Spitzer /IR AC, Spitzer /MIPS or WISE0 im- 
ages) and other observational evidence reported in the liter- 
ature (e.g., molecular outflows, masers, HII regions, radio- 
jets; see references in Table IBT]) , summarized in Table IBTI 
Within the starless cores we consider two groups: 'quiescent 



starless' and 'perturbed starless' cores. The former group 
includes cores appearing dark in the IRAC images, away 
(>10") from strong IR sources or outflows, with or with- 
out millimetre continuum. On the other hand, perturbed 
starless cores are neither associated with a point infrared 
source but lie close (<10") to a bright IR source and/or 
along the axis of a molecular outflow, suggesting that these 
cores might be perturbeqj. A core is classified as 'proto- 
stellar' if it is associated with an infrared point source, and 
with other star formation signposts (e.g., maser, outflow, 
centimetre continuum). Some cores were tentatively classi- 
fied as 'protostellar?' if the only hint of star formation is the 
association with an infrared point source, which could be a 
background or foreground star not associated with the dense 
core. These tentatively classified cores are not included in 
our statistical analysis when comparing the different types 
of cores, to avoid contamination by misclassined cores. In 
Fig. [1] the stars indicate the position of centimetre contin- 
uum sources, likely tracing HII regions or thermal radio jets. 
In many cases, the most massive source of the region, usu- 
ally associated with a strong centimetre source, tracing an 
(ultracompact) HII region is located offset of the core peaks 
(e.g., 20081+3122, 20293+3952, 20343+4129). 



2 The Wide-filed Infr ared Survey Explorer (WISE) is described 
in dWright et al.ll2010h . 



3 The differentiation between quiescent and pertur bed starless 
cores was already introduced in lFontani et al.l (|201ll) . 



Table 3. Physical parameters of the ammonia cores identified toward the 15 massive star-forming regions 



C^ 





# 


Offset a 

Ax Ay 


Size b 
(pc) 




NH 3 ( 


l,l) c 






NH 3 (2,2) c 




T rot e 


A^NH 3 e 


^ 


ID Region 


A x r 


V 


Av 


T* 


^mb 


V 


Av 


CO 


01 00117+6412 


1 


-20.5 


-7.0 


0.085 


8.2 + 1.0 


-36.2 + 0.1 


0.8 + 0.1 


1.80 + 0.29 


1.4 + 0.2 


-36.1 + 0.1 


1.0 + 0.2 


14.6 + 0.9 


8.4 + 0.9 


S3 


02 


2 


-12.5 


-16.5 


0.046 


4.5 + 1.3 


-35.5 + 0.1 


1.2 + 0.2 


1.40 + 0.62 


0.5 + 0.2 


-35.6 + 0.6 


2.4+1.3 


12.3 + 1.5 


10.1 + 2.3 


i 


03 


3 


-10.0 


-2.5 


0.039 


2.5 + 0.3 


-36.4 + 0.1 


1.6 + 0.3 


< 0.26 


0.6 + 0.2 


-36.0 + 0.5 


2.2 + 1.0 


15.8 + 0.8 


2.9 + 0.6 


04 


4 


-8.5 


-19.0 


0.026 


3.5 + 1.4 


-35.1 + 0.2 


1.4 + 0.3 


0.77 + 0.82 


0.8 + 0.3 


-35.0 + 0.3 


1.0 + 0.6 


15.4 + 2.9 


6.4 + 2.0 


i 


05 


5 


-23.0 


-14.5 


0.048 


8.1 + 2.3 


-36.3 + 0.1 


0.6 + 0.1 


3.45 + 0.97 


0.6 + 0.3 


-36.3 + 0.3 


0.8 + 0.6 


11.3 + 1.8 


13.2 + 3.4 


06 AFGL5142 


1 


+0.0 


-1.0 


0.040 


12.3 + 1.7 


-2.2 + 0.2 


3.8 + 0.2 


1.35 + 0.30 


7.1 + 0.6 


-4.0 + 0.2 


4.6 + 0.5 


36.0 + 3.9 


41.1 + 6.0 


53 

<0 


07 


2 


-6.0 


-3.5 


0.059 


12.9 + 1.5 


-3.2 + 0.1 


1.7 + 0.1 


0.69 + 0.25 


9.2 + 0.5 


-3.8 + 0.1 


1.4 + 0.1 


32.5 + 3.1 


16.5 + 2.3 


a> 


08 


3 


-12.0 


-19.0 


0.073 


9.2 + 2.2 


-5.0 + 0.1 


1.4 + 0.2 


1.19 + 0.53 


3.5 + 0.2 


-4.5 + 0.1 


1.6 + 0.1 


21.2 + 3.5 


12.4 + 2.7 


c-+- 


09 05345+3157NE 


1 


+6.4 


-8.8 


0.090 


11.8 + 0.6 


-17.7 + 0.1 


1.3 + 0.1 


1.64 + 0.13 


2.3 + 0.2 


-17.7 + 0.1 


1.7 + 0.1 


15.3 + 0.4 


16.4 + 0.9 


£- 


10 


2 


+0.0 


-4.0 


0.046 


10.8 + 0.9 


-18.3 + 0.1 


1.0 + 0.1 


1.06 + 0.17 


3.2 + 0.3 


-18.4 + 0.1 


0.9 + 0.1 


18.1 + 0.8 


9.0 + 0.8 




11 


3 


+2.4 


-21.6 


0.065 


17.7+1.1 


-18.1 + 0.1 


0.5 + 0.1 


1.90 + 0.16 


2.9 + 0.3 


-18.0 + 0.1 


0.5 + 0.1 


14.5 + 0.5 


10.2 + 0.7 




12 


4 


+5.6 


-33.6 


0.059 


8.8 + 0.8 


-17.1 + 0.1 


1.0 + 0.1 


1.27 + 0.21 


2.0 + 0.2 


-17.2 + 0.1 


1.1 + 0.2 


16.1 + 0.8 


10.8 + 0.8 




13 


5 


+8.0 


+4.0 


0.045 


17.3 + 1.4 


-17.6 + 0.1 


0.5 + 0.1 


1.95 + 0.19 


2.3 + 0.3 


-17.6 + 0.1 


0.8 + 0.1 


13.2 + 0.6 


10.2 + 0.9 




14 


6 


+11.2 


+9.6 


0.047 


17.4 + 1.3 


-17.9 + 0.1 


0.4 + 0.1 


2.53 + 0.23 


1.1 + 0.4 


-17.8 + 0.1 


0.6 + 0.2 


10.7 + 0.4 


12.0 + 1.0 




15 


7 


-20.8 


+2.4 


0.056 


10.9 + 1.1 


-17.9 + 0.1 


0.4 + 0.1 


1.48 + 0.23 


2.5 + 0.4 


-17.9 + 0.1 


0.4 + 0.1 


16.3 + 1.0 


5.1 + 0.5 




16 05358+3543NE 


1 


-3.2 


+6.0 


0.089 


9.5 + 0.7 


-34.6 + 0.1 


1.2 + 0.1 


< 0.14 


5.7 + 0.3 


-34.5 + 0.1 


1.4 + 0.1 


23.4 + 0.9 


7.1 + 0.8 




17 


2 


-14.0 


+8.0 


0.061 


14.5 + 1.5 


-35.4 + 0.1 


0.8 + 0.1 


1.14 + 0.22 


4.3 + 0.4 


-35.4 + 0.1 


1.3 + 0.1 


18.2 + 1.1 


9.9 + 1.1 




18 


3 


+2.0 


-2.8 


0.055 


5.9 + 0.3 


-35.0 + 0.1 


1.7 + 0.1 


< 0.07 


4.4 + 0.4 


-35.1 + 0.1 


1.8 + 0.2 


26.6 + 1.0 


6.2 + 0.6 




19 


4 


+4.8 


-18.8 


0.074 


6.7+1.2 


-35.8 + 0.1 


1.3 + 0.1 


0.30 + 0.36 


3.8 + 0.3 


-35.7 + 0.1 


1.2 + 0.1 


23.9 + 2.5 


6.5 + 1.2 




20 


5 


-24.4 


+9.6 


0.069 


16.5 + 3.0 


-36.1 + 0.1 


0.5 + 0.1 


2.26 + 0.50 


3.6 + 0.7 


-36.1 + 0.1 


0.7 + 0.2 


16.9 + 2.2 


9.5 + 1.7 




21 


6 


-9.6 


-23.6 


0.031 


15.6 + 2.6 


-36.7 + 0.1 


0.5 + 0.1 


1.69 + 0.38 


2.1 + 0.4 


-36.8 + 0.2 


1.3 + 0.3 


13.2 + 1.1 


8.7+1.4 




22 


7 


-21.2 


-45.2 


0.078 


7.6 + 1.3 


-37.3 + 0.1 


0.9 + 0.1 


1.10 + 0.36 


1.3 + 0.4 


-37.3 + 0.2 


1.1 + 0.4 


14.2 + 1.1 


12.9 + 1.2 




23 05373+2349 


1 


-9.6 


+8.8 


0.081 


12.3 + 0.5 


+2.0 + 0.1 


1.0 + 0.1 


1.23 + 0.13 


3.4 + 0.1 


+2.1 + 0.1 


1.4 + 0.1 


17.8 + 0.5 


10.4 + 0.6 




24 










18.3 + 1.6 


+3.2 + 0.1 


0.6 + 0.1 


3.15 + 0.28 


2.9 + 0.1 


+3.3 + 0.1 


0.7 + 0.1 


14.9 + 1.0 


13.2 + 1.3 




25 


2 


+0.0 


-4.8 


0.044 


11.1 + 0.6 


+1.5 + 0.1 


1.2 + 0.1 


1.16 + 0.16 


3.4 + 0.1 


+2.1 + 0.1 


1.4 + 0.1 


18.7 + 0.6 


11.2 + 0.6 




26 










9.7+1.5 


+2.6 + 0.1 


0.5 + 0.1 


0.57 + 0.30 


4.2 + 0.2 


+2.5 + 0.1 


0.5 + 0.1 


21.1 + 1.9 


3.4 + 0.6 




27 19035+0641 


1 


-0.5 


+0.3 


0.044 


5.9 + 1.1 


+32.4 + 0.2 


4.4 + 0.3 


0.61 + 0.36 


3.3 + 0.2 


+33.0 + 0.2 


5.0 + 0.4 


25.0 + 3.1 


21.0 + 3.8 




28 


2 


-5.8 


+0.3 


0.050 


4.7 + 0.3 


+32.9 + 0.1 


3.1 + 0.2 


< 0.15 


3.4 + 0.2 


+33.1 + 0.1 


3.0 + 0.2 


26.0 + 1.0 


9.2 + 0.9 




29 


3 


-4.5 


+11.5 


0.028 


5.1 + 0.9 


+34.2 + 0.1 


2.7 + 0.2 


0.30 + 0.33 


1.4 + 0.3 


+34.2 + 0.3 


2.3 + 0.4 


16.4 + 1.2 


11.7 + 2.0 




30 


4 


+9.5 


-0.3 


0.036 


4.0 + 0.8 


+32.2 + 0.2 


2.8 + 0.3 


0.33 + 0.43 


1.3 + 0.2 


+32.3 + 0.3 


2.8 + 0.5 


17.9 + 1.7 


9.1 + 1.8 




31 


5 


-15.5 


-6.0 


0.034 


6.4 + 1.3 


+31.9 + 0.1 


1.6 + 0.2 


1.52 + 0.46 


1.2 + 0.3 


+31.7 + 0.2 


1.3 + 0.4 


15.0 + 1.7 


13.2 + 2.3 




32 20081+3122 


1 


+1.6 


+15.2 


0.232 


15.6 + 1.7 


+10.5 + 0.1 


2.1 + 0.1 


2.10 + 0.28 


4.0 + 0.3 


+10.4 + 0.1 


2.7 + 0.2 


18.2 + 1.5 


31.1 + 3.4 




33 


2 


-16.8 


+9.6 


0.140 


7.6 + 0.9 


+11.9 + 0.1 


1.2 + 0.1 


0.89 + 0.36 


2.4 + 0.2 


+11.9 + 0.1 


1.4 + 0.2 


18.4 + 1.2 


8.7+1.0 




34 










6.4 + 0.9 


+13.6 + 0.1 


1.3 + 0.1 


0.23 + 0.30 


3.0 + 0.2 


+13.5 + 0.1 


1.6 + 0.2 


20.9 + 1.6 


6.5 + 1.0 




35 20126+4104 


1 


+0.0 


+0.8 


0.033 


21.9 + 2.4 


-4.3 + 0.2 


4.7 + 0.2 


2.09 + 0.28 


8.4 + 0.6 


-4.1 + 0.2 


5.5 + 0.5 


27.7 + 2.5 


81.5 + 9.6 




36 


2 


-4.8 


-12.0 


0.076 


85.7 + 5.0 


-4.5 + 0.1 


0.8 + 0.1 


5.40 ± 0.34 


11.8 + 0.4 


-4.5 + 0.1 


1.2 + 0.1 


15.7+1.0 


58.8 + 5.0 




37 


3 


-13.6 


-8.0 


0.084 


52.0 + 4.0 


-4.5 + 0.1 


0.9 + 0.1 


3.89 + 0.26 


9.3 + 0.3 


-4.5 + 0.1 


1.2 + 0.1 


17.0 + 1.3 


39.3 + 3.8 




38 G75.78+0.34 


1 


+0.7 


-0.7 


0.062 


5.1 + 0.6 


+0.2 + 0.3 


4.6 + 0.6 


< 0.22 


5.6 + 0.6 


+0.4 + 0.3 


4.5 + 0.6 


35.2 + 3.3 


15.7 + 2.8 




39 


2 


+17.5 


+7.7 


0.113 


5.2 + 0.4 


+0.4 + 0.2 


3.0 + 0.3 


< 0.11 


3.6 + 0.4 


+0.4 + 0.2 


2.8 + 0.3 


25.4 + 1.2 


10.5 + 1.2 




40 


3 


+18.9 


+11.9 


0.077 


6.3 + 0.8 


+1.1 + 0.2 


2.6 + 0.4 


< 0.27 


4.2 + 0.4 


+0.9 + 0.2 


2.8 + 0.3 


25.1 + 1.9 


11.1 + 1.9 




41 


4 


+23.8 


+13.3 


0.102 


5.6 + 0.3 


+1.5 + 0.1 


1.3 + 0.2 


0.22 + 0.03 


3.3 + 0.3 


+1.5 + 0.2 


1.7 + 0.3 


24.1 + 0.7 


5.8 + 0.6 




42 










5.7 + 0.7 


-0.2 + 0.1 


1.3 + 0.2 


< 0.20 


2.9 + 0.3 


-0.5 + 0.2 


1.8 + 0.3 


21.5 + 1.3 


5.4 + 0.9 




43 


5 


+27.3 


+14.7 


0.080 


6.3 + 0.9 


+0.3 + 0.1 


2.4 + 0.2 


0.14 + 0.37 


3.2 + 0.4 


+0.3 + 0.2 


3.0 + 0.4 


21.8 + 1.6 


12.4 + 1.7 




44 


6 


+35.7 


+11.9 


0.162 


8.6 + 1.2 


-0.2 + 0.1 


1.4 + 0.1 


0.17 + 0.36 


4.5 + 0.4 


-0.2 + 0.1 


1.3 + 0.1 


22.3 + 1.8 


10.8 + 1.4 





Table 3 - continued 





# 


Offset a 

Ax Ay 


Size b 
(pc) 




NH 3 ( 


l,l) c 






NH 3 (2,2) c 




T ro t e 




ID Region 


Ax T 


V 


Av 


r d 


^mb 


V 


Av 


A^NH 3 e 


45 20293+3952 


1 


+0.0 


+18.0 


0.089 


16.7 + 4.4 


+7.3 + 0.1 


1.1 + 0.1 


3.33 + 0.88 


3.0 + 0.6 


+7.4 + 0.2 


1.3 + 0.3 


16.3 + 3.6 


20.4 + 5.5 


46 


2 


+36.0 


-10.0 


0.148 


23.1 + 3.9 


+6.0 + 0.1 


0.9 + 0.1 


4.84 + 0.72 


2.6 + 0.3 


+6.1 + 0.1 


1.6 + 0.2 


13.4 + 2.0 


35.4 + 5.8 


47 


3 


+32.0 


+8.0 


0.112 


6.4 + 1.0 


+5.9 + 0.1 


2.1 + 0.2 


1.63 + 0.37 


2.0 + 0.2 


+6.0 + 0.2 


2.9 + 0.4 


20.1 + 2.2 


18.7 + 2.1 


48 


4 


+29.0 


-2.5 


0.125 


13.3 + 2.3 


+5.5 + 0.1 


1.0 + 0.1 


2.87 + 0.55 


2.4 + 0.3 


+5.6 + 0.1 


1.8 + 0.2 


15.6 + 2.0 


18.5 + 2.8 


49 


5 


+15.5 


-0.5 


0.083 


12.8 + 2.8 


+5.3 + 0.1 


1.0 + 0.1 


2.80 + 0.63 


2.7 + 0.3 


+5.4 + 0.1 


1.6 + 0.2 


17.5 + 3.2 


15.0 + 3.2 


50 


6 


-24.5 


+40.5 


0.103 


8.5 + 1.8 


+4.9 + 0.1 


1.3 + 0.1 


2.18 + 0.54 


1.4 + 0.2 


+4.9 + 0.2 


1.9 + 0.3 


14.8 + 1.9 


23.9 + 2.8 


51 20343+4129 


1 


+12.9 


+0.8 


0.044 


11.6 + 2.3 


+11.8 + 0.1 


1.9 + 0.2 


2.04 + 0.52 


2.6 + 0.6 


+11.9 + 0.3 


2.5 + 0.6 


16.7 + 2.3 


23.9 + 4.6 


52 


2 


-9.1 


-0.8 


0.035 


6.7 + 0.9 


+10.9 + 0.1 


1.6 + 0.3 


< 0.24 


2.5 + 0.7 


+10.9 + 0.2 


0.9 + 0.3 


18.5 + 1.2 


7.1 + 1.6 


53 


3 


+6.8 


+8.0 


0.023 


4.6 + 1.1 


+13.0 + 0.2 


1.3 + 0.4 


< 0.25 


2.2 + 1.7 


+12.4 + 0.4 


0.4 + 0.3 


20.8 + 2.5 


3.8 + 1.5 


54 


4 


+1.1 


-0.4 


0.036 


6.1 + 2.0 


+10.2 + 0.1 


1.1 + 0.2 


0.60 + 0.67 


2.1 + 0.5 


+10.2 + 0.2 


1.0 + 0.3 


18.8 + 3.3 


5.3 + 1.8 


55 22134+5834 


1 


-10.4 


-20.7 


0.071 


9.7+1.8 


-19.2 + 0.1 


0.7 + 0.1 


3.15 + 0.57 


1.2 + 0.2 


-19.2 + 0.1 


0.8 + 0.2 


13.1 + 1.6 


12.2 + 2.1 


56 


2 


+18.6 


-11.0 


0.094 


5.1 + 0.5 


-18.7 + 0.1 


1.0 + 0.1 


1.21 + 0.22 


2.1 + 0.2 


-18.7 + 0.1 


1.0 + 0.1 


22.7+1.7 


5.7 + 0.5 


57 


3 


-17.9 


-11.0 


0.107 


7.4 + 0.6 


-18.7 + 0.1 


0.6 + 0.1 


1.94 + 0.22 


1.1 + 0.2 


-18.7 + 0.1 


0.8 + 0.2 


13.9 + 0.7 


7.1 + 0.6 


58 


4 


+22.8 


-14.5 


0.100 


3.2 + 0.6 


-18.4 + 0.1 


1.2 + 0.1 


0.49 + 0.35 


1.0 + 0.1 


-18.4 + 0.1 


1.2 + 0.2 


17.3 + 1.6 


4.1 + 0.6 


59 


5 


+9.0 


-17.3 


0.068 


3.5 + 0.7 


-18.7 + 0.1 


0.9 + 0.1 


0.92 + 0.41 


1.2 + 0.2 


-18.9 + 0.1 


1.0 + 0.2 


19.7 + 2.3 


3.6 + 0.5 


60 22172+5549N 


1 


+4.5 


+0.3 


0.099 


4.2 + 0.6 


-43.2 + 0.1 


1.4 + 0.1 


0.25 + 0.26 


1.8 + 0.2 


-43.3 + 0.1 


1.5 + 0.2 


20.3 + 1.4 


4.5 + 0.6 


61 


2 


+3.6 


+9.9 


0.040 


5.6 + 1.5 


-43.9 + 0.1 


0.9 + 0.1 


1.06 + 0.57 


1.3 + 1.7 


-44.0 + 2.0 


0.8 + 3.0 


16.0 + 2.2 


5.3 + 1.3 


62 


3 


+4.2 


-6.0 


0.039 


3.2 + 0.8 


-43.6 + 0.1 


1.8 + 0.2 


1.31 + 0.54 


0.9 + 0.6 


-43.5 + 1.0 


2.3 + 1.0 


17.5 + 2.6 


8.2 + 1.5 


63 22198+6336 


1 


-11.0 


+11.0 


0.023 


11.2 + 2.2 


-10.6 + 0.1 


1.0 + 0.1 


1.38 + 0.43 


1.5 + 0.6 


-10.7 + 0.2 


1.0 + 0.5 


13.2 + 1.2 


12.9 + 2.5 


64 


2 


-14.5 


+24.3 


0.021 


7.2 + 1.9 


-10.8 + 0.1 


0.9 + 0.1 


0.83 + 0.50 


1.9 + 0.9 


-10.8 + 0.1 


0.5 + 0.2 


16.6 + 2.1 


7.3 + 1.6 


65 


3 


+0.3 


+4.5 


0.021 


6.4 + 1.4 


-10.5 + 0.1 


1.1 + 0.1 


0.60 + 0.40 


2.4 + 0.4 


-10.7 + 0.2 


1.4 + 0.3 


19.4 + 2.2 


5.9 + 1.2 


66 


4 


-5.3 


+2.3 


0.031 


9.2 + 1.4 


-10.6 + 0.1 


0.9 + 0.1 


1.56 + 0.33 


1.5 + 0.3 


-10.7 + 0.2 


1.9 + 0.5 


14.1 + 1.1 


10.2 + 1.4 


67 


5 


-17.0 


+17.3 


0.020 


16.1 + 3.6 


-10.8 + 0.1 


0.5 + 0.1 


2.33 + 0.62 


2.4 + 1.2 


-10.7 + 0.2 


0.4 + 0.2 


14.1 + 1.9 


9.0 + 2.1 


68 


6 


-9.8 


-4.5 


0.017 


3.8 + 0.5 


-10.9 + 0.1 


1.2 + 0.2 


< 0.29 


1.6 + 0.6 


-11.1 + 0.2 


0.9 + 0.5 


19.4 + 1.2 


3.2 + 0.7 


69 


7 


+6.3 


-9.0 


0.024 


11.5 + 2.9 


-11.5 + 0.1 


0.6 + 0.1 


2.52 + 0.70 


2.5 + 0.6 


-11.4 + 0.1 


0.8 + 0.3 


17.1 + 3.3 


8.2 + 2.1 


70 


8 


-27.3 


+32.5 


0.023 


16.3 + 3.9 


-11.5 + 0.1 


0.3 + 0.1 


2.93 + 0.73 












71 CepA 


1 


+41.6 


+14.4 


0.039 


8.9 + 1.0 


-11.7 + 0.1 


1.3 + 0.1 


0.54 + 0.24 


4.1 + 0.3 


-11.6 + 0.1 


1.5 + 0.1 


21.9 + 1.5 


13.4 + 1.1 


72 


2 


+4.8 


+32.8 


0.041 


22.4 + 2.1 


-10.2 + 0.1 


0.5 + 0.1 


2.13 + 0.26 


5.7 + 0.5 


-10.2 + 0.1 


0.8 + 0.1 


18.3 + 1.3 


11.5 + 1.0 


73 


3 


+49.6 


+21.6 


0.027 


7.4 + 1.4 


-11.9 + 0.1 


1.5 + 0.1 


0.39 + 0.32 


3.1 + 0.3 


-12.0 + 0.1 


1.9 + 0.2 


20.2 + 2.0 


15.3 + 1.6 


74 


4 


+32.0 


+6.4 


0.036 


9.7 + 0.8 


-12.0 + 0.1 


0.9 + 0.1 


0.48 + 0.18 


5.2 + 0.3 


-12.0 + 0.1 


1.0 + 0.1 


23.8 + 1.3 


7.7 + 0.6 


75 


5 


-11.2 


-15.2 


0.056 


10.4 + 1.6 


-9.5 + 0.1 


0.8 + 0.1 


0.49 + 0.27 


6.3 + 0.3 


-9.5 + 0.1 


0.9 + 0.1 


26.0 + 2.8 


6.1 + 1.0 


76 


6 


+65.6 


+30.4 


0.041 


13.7+1.5 


-12.1 + 0.1 


0.7 + 0.1 


2.17 + 0.30 


3.2 + 0.2 


-12.1 + 0.1 


1.0 + 0.1 


17.5 + 1.4 


31.1 + 1.2 


77 


7 


+10.4 


+17.6 


0.023 


17.8 + 2.4 


-10.8 + 0.1 


0.5 + 0.1 


2.33 + 0.34 


5.1 + 0.6 


-10.8 + 0.1 


0.7 + 0.1 


20.6 + 2.3 


8.0+1.1 



a Offset (in arcsec) between the position of the core and the phase centre (see Table [2]). x and y correspond to right ascension and declination, respectively. 

b Linear diametre of the clump computed from the deconvolved angular diametre, 0g, which has been calculated assuming the cores are Gaussians from the expression 

Q2 _ t?\\tu\/\2 _ WP"RW2 



FWHM 2 — HPBW 2 , with FWHM = 2^A/tt and A being the area inside the polygons (shown in Fig. [T} and with HPBW corresponding to the geometric mean of the minor and 
major axes of the synthesised beams (see Table [2]). 
c Parameters obtained from the spectral fit to the NH3 (1,1) and (2,2) lines: v and Ai; are the velocity of the local standard of rest and the linewidth in km s — 1 , and T m b is the 
temperature of the line in K. The hyperfine structure of the (1,1) transition allows us to measure the opacity r of the main line. The parameter Axris defined as 
A = f [J u (T ex ) - J„(T hg )] with / = 1 (see Appendix[A}. 

d For a few cases, corresponding to the cores with optically thin emission, we provide an upper limit in the opacity. Even if r cannot be determined with high precision, the derived 
quantities, T ro t and A^nh 3 , m the optically thin limit, depend on the line intensity, A x r, which is very well constrained by the data (see Appendix [A}. 
e Rotational temperature (in K) and column density (in 10 14 cm -2 ). 
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3.2 Averaged physical parameters for ammonia 
dense cores 

For each core, we extracted the spectra averaged over the 
area assigned by the polygons (shown in Fig. [T]) in both 
NH 3 (1,1) and (2,2) lines, and fitted the hyperfine struc- 
ture of the NH3 (1,1) transition following the procedure de- 
scribed in AppendixlAl and a Gaussian to the NH3 (2,2) line 
(see observed and fitted spectra in Fig. IB2|) . This allowed 
us to obtain the opacity, velocity, and linewidth for each 
corqj. We derived the rotational temperature and primary 
beam corrected ammonia column density following the pro- 
cedure outlined in the appendix of iBusquet et al.l ( 20091 ). All 
these parameters are listed in Table [3] The average values of 
linewidth, opacity of the NH3 (1,1) main line, rotational tem- 
perature, and ammonia column density of all the 73 cores 
are ~ 1.4 km s _1 , 1.4, 19 K, and 10 15 cm -2 , respectively. 
Our values of linewidth, rotational temperature and ammo- 
nia cojmmi_dcnsity are v ery similar to the values obtained 
by lUrauhart et"atl (|201lh for the RMS (Red MSX Source) 
survey, of 2 km s -1 , 22 K and 2 x 10 15 cm -2 , even using 
very different observing techniques and sampling different 
spatial scales (single-dish vs interferometry) . If we analyse 
the distributions of linewidth, opacity, rotational tempera- 
ture, and column density for the starless and protostellar 
populations separately (see Fig. [2] panels A to E), the aver- 
age values that we obtain are (i) ~ 1.0 km s _1 , 1.5, 16.0 K, 
and 10 15 cm -2 for quiescent starless cores, (ii) ~ 1.4 km s _1 , 
1.2, 19.3 K, and 10 15 cm -2 for perturbed starless cores, and 
(hi) - 1.8 km s"\ 1.2, 21.3 K, and 10 15 cm" 2 for protostel- 
lar cores. Thus, quiescent starless cores have, as expected, 
linewidths and rotational temperatures smaller than pro- 
tostellar cores (see Fig. [2] panels A and D; Table 0. The 
linewidths and temperatures of perturbed starless cores are 
in between the values measured for quiescent starless and 
protostellar cores, consistent with the fact that these cores 
were originally quiescent and then perturbed by the pas- 
sage of outflows or the high radiation of a nearby bright 
infrared source. These co res might be simila r to the 'warm' 
starless cores reported in lMotte et al.l (|2010l ) using Herschel 
observations towards Cygnus X, and f or some of them a low 
deuterium fractionation is measured (jFontani et al.ll201lf ). 
Further studies should be carried out to assess the true na- 
ture of these perturbed starless cores. 

The average values derived for the rotational tempera- 
ture of quiescent starless and protostellar cores (of 16.0 K 
and 21.3 K, respectively, corresponding to kinetic temper - 
atures of 18.8 K and 28.8 K, following Pfafalla et al.ll2004l ) 
can be compared to the values reported in the literature, 
although derived f r om ob servations of lower angular reso- 
lution. lOlmi et al.l (|2013T ) studied the dense clumps in the 
Galactic field / = 59° (at a distance of 3.6 kpc, similar to 
the dista nces of our regions) within the Hi-GAL Herschel 
project (jMolinari et al.ll2010l h and derived kinetic temper- 
atures of 17 K for starless cores and 20 K for protostellar 
cores, following the same trend as we found in our cores. 
Urquhart et al. (2011) estimated the kinetic temperatures 
for clumps associated with young stellar objects and HII re- 
gions in the RMS survey to be in the range 20-25 K, also in 



4 For four cores it was necessary to fit the ammonia spectra with 
two different velocity components (see Table [3] and Fig. IB2|) . 



agreement with our protostellar temperatures. However, it 
is worth noting that the average rotational temperature of 
16 K found in the quiescent starless cores of our clustered 
massive star-forming regions is significantly larger than the 
values measured in starless cores o f low-mass isolated regions 
(of about 10 K or even less; e. g.. iTafalla et alll200l [20o3; 
JRathborne et al.ll2008l : iRosolowskv et al.ll2008l : ISchnee et"al] 
2010), and of infrared dark clouds with no active star for- 
mation (rsj 10 K; e.g., Rag an et al.ll2011l ) . Temperatures of 
~ 15 K for starless cores are also found by I Li et al] (|2003r ) in 
the quiescent cores near the Trapezium cluster. Li et al. find 
that the temperatures of the massive quiescent cores can be 
well explained by the dust being heated by the external UV 
field from the Trapezium cluster, at a distance of ~ 1 pc 
from the massive cores. This possibility is investigated for 
the regions in our sample in Sect. 14. ll 

The linewidths measured in our starless cores have an 
important non-thermal component, whic h is almost negligi- 
ble in low-ma ss clustered regions fe. g. jFriesen et al.l l2009; 
iFoster et al.l l2QQ9h . To evaluate this component, we as- 
sumed the relation cr n th = V^obs — a th wnere ^obs is the 
measured velocity dispersion of the core (= Av/(81n2) 1//2 
for a Gaussian line profile with Av the measured FWHM 
of the line) and <7th is the thermal velocity dispersion. For 
NH3, assuming a Maxwellian velocity distribution, cr t h can 
be calculated from y/k^ Tki n /(/i ^h), with /cb the Boltz- 
mann constant, \i the molecular weight (= 17.03 for am- 
monia), ??iH the mass of the Hydrogen atom, and T^ n the 
kinetic temperature listed in Col. 3 of Table [H In Cols. 5 
and 6 of the same Table we list the non-thermal velocity dis- 
persion (cr n th) and the non-thermal velocity dispersion over 
the isothermal sound speed (cr n th/c s ), where c s is calculated 
using a mean molecular weight \i — 2.33. From the distri- 
butions shown in Fig. [2] (panels H and I) we see that pro- 
tostellar cores have a non-thermal component (0.76 km s -1 , 
or equivalently Av n th = 1.8 km s _1 ) higher than the non- 
thermal component of the quiescent starless (0.44 km s -1 , 
or A^nth = 1.0 kms) and perturbed starless (0.56 km s _1 , 
or A^nth = 1.3 km s _1 ) cores. For the cr n th/c s we obtain 
mean ratios of 1.7, 1.9 and 2.2 for quiescent starless, per- 
turbed starless and protostellar cores, respectively (see Ta- 
ble [5]), indicating that the cores are characterised by moder- 
ately supersonic non- thermal motions. In general our cores 
(both starless and protostellar), associated with intermedi- 
ate and high-mass star-forming regions, have values of tur- 
bulent Mach number ~ 2, being greater than those found in 
Perseus or Ophiucus (with values between 0.5 and 1.5; e . g . , 
I Andre et aJ1l2007l : iFriesen et alll2009l : IFoster et al]l20Q9h . 

From the ammonia column density (listed in Table [3}, 
we derived the mass of the core and the core mass surface 
density, following the procedures explained in the footnotes 
of Table 2) These two quantities are listed in Cols. 8 and 
9 of Table [H and their distributions are shown in panels F 
and G of Fig. [2] Starless and protostellar cores show similar 
mean masses (~ 1.0-1.8 Mq) and mean surface densities 
(~ 0.10-0.12 g cm" 2 ) (see Table©. 



3.3 Ammonia correlations 

As presented in Sect. 13.21 we note some differences in the pa- 
rameters derived for the starless and protostellar cores. In 
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Table 4. Derived parameters for the ammonia cores 







7kin b 


^th 


^nth 




Ctot 


M d 

i\a core 


v e 

^core 


M vir f 


ID 


Type a 


(K) 


(km s- 1 ) 


(km s- 1 ) 


^nth/Cs 


(kms- 1 ) 


(M ) 


(g cm -2 ) 


(M Q ) 


01 


P 


16.4 


0.09 


0.32 


1.32 


0.40 


2.56 


0.09 


5.42 


02 


S 


13.3 


0.08 


0.51 


2.36 


0.56 


0.90 


0.11 


7.19 


03 


P 


18.1 


0.09 


0.67 


2.63 


0.71 


0.18 


0.03 


10.26 


04 


s* 


17.5 


0.09 


0.60 


2.41 


0.65 


0.18 


0.07 


5.59 


05 


p? 


12.1 


0.08 


0.26 


1.25 


0.33 


1.28 


0.15 


2.03 


06 


p 


66.3 


0.18 


1.59 


3.27 


1.66 


2.76 


0.46 


59.31 


07 


p? 


53.9 


0.16 


0.71 


1.61 


0.83 


2.41 


0.18 


17.99 


08 


p? 


26.9 


0.11 


0.60 


1.94 


0.67 


2.78 


0.14 


15.76 


09 


p 


17.4 


0.09 


0.56 


2.27 


0.62 


5.57 


0.18 


17.10 


10 


p 


21.6 


0.10 


0.41 


1.47 


0.49 


0.80 


0.10 


4.69 


11 


s 


16.2 


0.09 


0.21 


0.87 


0.32 


1.81 


0.11 


1.93 


12 


s 


18.5 


0.09 


0.40 


1.56 


0.47 


1.58 


0.12 


5.80 


13 


s 


14.5 


0.08 


0.21 


0.92 


0.31 


0.86 


0.11 


1.32 


14 


s 


11.3 


0.07 


0.17 


0.83 


0.26 


1.11 


0.13 


0.90 


15 


s 


18.9 


0.10 


0.16 


0.63 


0.31 


0.67 


0.06 


1.16 


16 


p 


31.1 


0.12 


0.49 


1.48 


0.59 


2.36 


0.08 


13.23 


17 


p? 


21.8 


0.10 


0.33 


1.18 


0.43 


1.54 


0.11 


4.20 


18 


p 


38.0 


0.14 


0.69 


1.89 


0.78 


0.79 


0.07 


16.01 


19 


s* 


32.1 


0.12 


0.54 


1.60 


0.64 


1.50 


0.07 


13.19 


20 


p 


19.7 


0.10 


0.20 


0.77 


0.33 


1.90 


0.11 


2.05 


21 


s* 


14.6 


0.08 


0.17 


0.76 


0.28 


0.35 


0.10 


0.66 


22 


s 


15.9 


0.09 


0.37 


1.55 


0.44 


3.28 


0.14 


6.49 


23 


p 


21.1 


0.10 


0.40 


1.45 


0.48 


2.87 


0.12 


7.92 


24 


p 


16.8 


0.09 


0.25 


1.04 


0.35 


3.62 


0.15 


3.42 


25 


s* 


22.5 


0.10 


0.49 


1.75 


0.57 


0.91 


0.13 


6.54 


26 


s* 


26.6 


0.11 


0.17 


0.56 


0.35 


0.28 


0.04 


1.10 


27 


p 


34.4 


0.13 


1.87 


5.37 


1.91 


1.70 


0.23 


90.38 


28 


s* 


36.6 


0.13 


1.30 


3.62 


1.35 


0.97 


0.10 


49.97 


29 


s 


18.9 


0.10 


1.16 


4.48 


1.19 


0.38 


0.13 


22.09 


30 


s* 


21.3 


0.10 


1.17 


4.24 


1.20 


0.49 


0.10 


28.71 


31 


s 


16.9 


0.09 


0.67 


2.72 


0.71 


0.64 


0.15 


8.92 


32 


p? 


21.8 


0.10 


0.89 


3.21 


0.93 


70.21 


0.35 


108.66 


33 


s 


22.0 


0.10 


0.51 


1.83 


0.58 


7.19 


0.10 


22.24 


34 


s 


26.2 


0.11 


0.56 


1.84 


0.64 


5.31 


0.07 


26.71 


35 


p 


40.5 


0.14 


1.98 


5.23 


2.01 


3.72 


0.91 


75.60 


36 


s 


17.9 


0.09 


0.32 


1.28 


0.41 


14.24 


0.66 


4.96 


37 


p? 


19.8 


0.10 


0.36 


1.35 


0.45 


11.64 


0.44 


6.77 


38 


p 


63.2 


0.18 


1.93 


4.07 


1.98 


2.54 


0.18 


135.01 


39 


s* 


35.3 


0.13 


1.28 


3.62 


1.33 


5.61 


0.12 


109.15 


40 


p? 


34.5 


0.13 


1.09 


3.12 


1.15 


2.76 


0.12 


54.11 


41 


s 


32.5 


0.13 


0.55 


1.61 


0.64 


2.55 


0.07 


18.66 


42 


s 


27.4 


0.12 


0.54 


1.72 


0.62 


2.35 


0.06 


17.82 


43 


s 


27.9 


0.12 


1.03 


3.26 


1.07 


3.33 


0.14 


49.72 


44 


p? 


29.0 


0.12 


0.58 


1.82 


0.66 


11.88 


0.12 


33.29 


45 


s* 


18.9 


0.10 


0.44 


1.71 


0.51 


6.78 


0.23 


10.54 


46 


s* 


14.8 


0.08 


0.36 


1.57 


0.43 


32.56 


0.40 


11.82 


47 


p 


24.9 


0.11 


0.90 


3.02 


0.94 


9.81 


0.21 


53.15 


48 


p? 


17.8 


0.09 


0.42 


1.68 


0.49 


12.15 


0.21 


13.55 


49 


p 


20.6 


0.10 


0.43 


1.59 


0.51 


4.33 


0.17 


9.43 


50 


p? 


16.6 


0.09 


0.56 


2.29 


0.61 


10.63 


0.27 


18.96 


51 


s* 


19.4 


0.10 


0.82 


3.14 


0.86 


1.94 


0.27 


17.51 


52 


s* 


22.1 


0.10 


0.66 


2.37 


0.72 


0.36 


0.08 


9.21 


53 


p? 


26.1 


0.11 


0.54 


1.77 


0.62 


0.08 


0.04 


4.03 


54 


p 


22.7 


0.10 


0.43 


1.53 


0.52 


0.29 


0.06 


4.18 


55 


s 


14.3 


0.08 


0.29 


1.28 


0.37 


2.59 


0.14 


3.75 


56 


s* 


29.7 


0.12 


0.40 


1.24 


0.52 


2.11 


0.06 


9.62 


57 


s 


15.4 


0.09 


0.26 


1.11 


0.35 


3.41 


0.08 


4.67 


58 


s 


20.4 


0.10 


0.48 


1.79 


0.55 


1.73 


0.05 


14.06 


59 


p? 


24.3 


0.11 


0.35 


1.21 


0.46 


0.70 


0.04 


5.42 
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Table 4 - continued 



ID Type a 



^kin 

(K) 



(km s- 1 ) 



c^nth 
(km s- 1 ) 



0~nth/Cs 



O~tot 

(km s- 1 ) 



M CI 



A A, 



(M ) (g cm- 2 ) (M ) 



60 


P 


25.2 


0.11 


0.60 


2.00 


0.67 


1.85 


0.05 


21.32 


61 


S 


18.3 


0.09 


0.36 


1.39 


0.44 


0.36 


0.06 


3.15 


62 


P? 


20.6 


0.10 


0.74 


2.75 


0.79 


0.52 


0.09 


12.75 


63 


s 


14.5 


0.08 


0.41 


1.80 


0.47 


0.29 


0.14 


2.32 


64 


p? 


19.3 


0.10 


0.39 


1.49 


0.47 


0.14 


0.08 


1.97 


65 


p 


23.7 


0.11 


0.46 


1.60 


0.55 


0.11 


0.07 


2.76 


66 


s* 


15.7 


0.09 


0.39 


1.67 


0.46 


0.41 


0.11 


2.92 


67 


s 


15.7 


0.09 


0.17 


0.73 


0.29 


0.15 


0.10 


0.43 


68 


s 


23.6 


0.11 


0.52 


1.79 


0.59 


0.04 


0.04 


2.78 


69 


s 


20.1 


0.10 


0.26 


0.96 


0.37 


0.20 


0.09 


1.05 


71 


s* 


28.2 


0.12 


0.55 


1.75 


0.64 


0.85 


0.15 


7.25 


72 


s* 


22.0 


0.10 


0.18 


0.65 


0.33 


0.81 


0.13 


1.04 


73 


s* 


25.1 


0.11 


0.61 


2.06 


0.68 


0.47 


0.17 


6.11 


74 


p? 


31.9 


0.12 


0.36 


1.07 


0.49 


0.42 


0.09 


3.05 


75 


s* 


36.6 


0.13 


0.30 


0.85 


0.47 


0.81 


0.07 


3.60 


76 


s* 


20.7 


0.10 


0.29 


1.07 


0.40 


2.19 


0.35 


2.25 


77 


p 


25.7 


0.11 


0.17 


0.57 


0.35 


0.18 


0.09 


0.56 



a Type of the core: S: quiescent starless, S*: perturbed starless, P: protostellar, and P?: possible protostellar (see Sect. I3.l1 and 

Table [BlJ. 

b Kinetic temperature, in K, derived from the rotational temperature and following Tafal la et al.l (|2004). We emphasise that Tafalla et 

al. empirical relation is accurate only for temperatures < 20 K, and thus the kinetic temperature given here should be regarded, for 

high temperatures, as an approximation. 

c Thermal (cr th ) and non-thermal (cr nt h) velocity dispersion, in km s _1 . 

d Mass of the core derived as M core = (N/X)/j, mnA, with /i (= 2.8) the mean molecular mass per H molecule, A (= 7r ( size/2) 2 ) the 

area, and assuming an abundance of NH3 with respect to H 2 (X) for all the cores of 4.2 x 10 -8 (e. g. JPillai et al.ll2006l : iFoster et al.l 

12009: iFriesen et al1l2009l: iRvgl et al.lkoiol: IChira et al.ll2013h . 

e Core mass surface density derived as S cor e = M core /A. 

f Virial mass derived from M v i r 



= 210 [R/pc] [At;/km s x ] 2 Mq assuming uniform density across the core (MacLare n et al.lll988l : 
lEstalella et al.l n.993). R is the radius of the core (obtained from the diametre size listed in Table [3j, and Av is the intrinsic linewidth of 
the NH 3 (1,1) line (listed in Table [3]). 



Table 5. Mean values derived from the NH3 observations (from Tables [3] and 0} 





T r ot 


Tkin 


size 


T 


^NH 3 


IVlcore 


^core 


A ^(l,l) 


Cnth 


Cnth /c s 


M vir 


Core type 


(K) 


(K) 


(pc) 




(cm" 2 ) 


(M ) 


(g cm~ 2 ) 


(km s- 1 ) 


(kms- 1 ) 




(M©) 


Quiescent starless 


16.0 


18.8 


0.063 


1.5 


10 15.0 


1.15 


0.10 


1.03 


0.44 


1.67 


5.1 


Perturbed starless 


19.3 


24.2 


0.056 


1.2 


10 15.0 


1.09 


0.12 


1.36 


0.56 


1.93 


7.0 


Protostellar 


21.3 


28.8 


0.063 


1.2 


10 15.0 


1.47 


0.12 


1.81 


0.76 


2.24 


11.7 


All sample 


18.8 


23.8 


0.065 


1.4 


10 15.0 


1.36 


0.12 


1.36 


0.57 


1.91 


7.8 



particular, an increase of the temperature and linewidth as 
we move toward more evolved objets (protostellar cores). In 
a search for a physical relation between the derived param- 
eters, we plotted several quantities against others, shown as 
six scatter plots in Fig. [3] In these plots we show the distri- 
bution of starless and protostellar core samples in red-circles 
and blue- stars, respectively. We performed some statistical 
tests to search for possible correlations between the physi- 
cal quantities. In the top right corner of each plot, we give 
Kendall's r rank correlation coefficient and the Pearson cor- 
relation coefficient (p) to a least-squares fit (dotted lines). 
We find that in five of the six plots there is a possible corre- 
lation, with r coefficients ~ 0.3, and p in the range 0.4-0.6. 

In the top panels of Fig.[3j we show the relation between 
the temperature, the linewidth and the mass of the core with 
the bolometric luminosity of the region, while in the bottom 



panels, we compare the temperature, the ammonia column 
density and the size of the core with the observed linewidth. 
We note that we are not able to derive a bolometric luminos- 
ity for each individual ammonia core but a single luminosity 
of the whole regiorjj. Thus, the plots in Fig. [3] must be inter- 
preted as a correlation between the properties of the cores 
within a region with respect to the luminosity of the re- 
gion (which is in general dominated by the brightest mid-IR 
source) . We find a correlation in five of the six panels (panels 
A to E of Fig. [3]). Similar trends between the same parame- 



5 To estimate the bolometric luminosity Lbol (listed in Ta- 
ble^), we compiled the flux densities at different frequencies from 
2MASS, Spitzer /IRAC+MIPS, WI SE, MSX, IRAS and JCMT 
(SCUBA; |Pi Francesco et IdJ |2008 [) ca talogues and followed the 
methodology explained in lPalau et al.l (J2013T ). 
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A 


. ^^ Quiesc. starless: 


24 


0.93 


1.03 


± 0.59 


15 


- 


~^~-~-j Pertur. starless: 


19 


1.25 


1.36 


± 0.81 






Protostellar: 


19 


1.27 


1.81 


±1.42 
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" E 


. ^^ Quiesc. starless: 24 

— — J Pertur. starless: 19 

Protostellar: 19 


1.4 1.5 ±1.3 
1 1.2 ± 1.3 
1.3 1.2 ± 1 ~ 
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L J 












Quiesc. starles 
Pertur. starles 
Protostellar: 



0.057 0.063 ± 0.037 
0.044 0.056 ± 0.033 
0.066 0.063 ± 0.027 



tzzzF^ 



t (opacity of main NH 3 (1,1) transition) 



0.05 0.1 0.15 0.2 

deconvolved diameter (pc) 



15 


- C 


• ^^ 


i 

Quiesc. starless: 
Pertur. starless: 
Protostellar: 


19 
19 


15.8 
18.6 
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16 

19.3 

21.3 


± 3.6 

± 4.1 " 
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Quiesc. starle: 
Pertur. starles 
Protostellar: 



rotational temperature (K) 



14.5 15 15. 5_ 

Log[ NH 3 column density (cm ) ] 





F. 
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Figure 2. Distributions of: (a) NH3 (1,1) linewidth, (b) opacity, (c) deconvolved diametre, (d) rotational temperature, (e) ammonia 
column density, (f) core mass M CO re, (g) core surface density Score, (h) non-thermal dispersion cr nth , and (i) non-thermal to thermal 
velocity dispersion ratio <r nth /c s , for quiescent starless (solid line histrogram), perturbed starless (dashed line histogram), and protostellar 
(grey filled histogram) cores. The numbers at the top of each panel are the total number of sources used in the histogram, the median, 
and the mean ± standard deviation values. 



ters have been reported from single-dish observatio nal stud- 
ie s fe.g.JChurchwell et"aIlll99Ql:lMvers et al.lll99ll: Jiiina et 



al.Fl 
all 



1999; Rvgl et al. 2010; Urauhart et al.ll201ll: Sepulveda et 
l201lh . Our correlations found in ammonia at higher an- 
gular resolution thus confirm that these relations take place 
also at the smaller spatial scales (~ 0.05 pc) of individual 
cores, and are not an observational (average) artifact due to 
the poor angular resolution of single-dish observations when 
studying clumps of sizes ~ 0.1-10 pc. It is worth noting that 
in panel B (linewidth vs luminosity) there are four cores 
with linewidths ~ 4 km s _1 , located well above the trend 
found for the other cores. These four cores (AFGL 5142-1; 
19035+0641-1; 20126+4104-1; G75.78+0.34-1) are known 
to be associated with powerful outflows that could perturb 
the motions of the dense gas close to the powering source 
fseelZhang et al.ll2002l: ICesaroni et alJ^OOdTSanchez-Mongi 
l201ll ; ISanchez-Monge et al.1 l2013l i~Finallv. in panel F, we 
investigate a possible relation between the size of the core 
and its linewidth. Since the publication of the Larson's scal- 
in g^ relations ( Larson 1981), several authors (e.g., M yers et 
al. Il983l ; ICaselli fe Mverdll995l ; IPirogov et alJl2003h found, 
from single-dish studies, that larger dense cores have broader 
lines, as expected in virialised cores. In our sample, however, 
we find no correlation (with r — —0.03 and p — 0.01) be- 
tween these two quantities. A possible explanation for this 



could be that a number of cores suffer an important injection 
of turbulence, e.g., from the passage of powerful outflows, 
which should affect the gas motions of t he cores only at the 
scales of our sample (^ 0.05 pc; see e. g. , IPalau et al.ll2007bl ; 
iLiu et al.|[2010l : iNakamura fe Lil l201lh . and is diluted at 
the larger scales studied in single-dish observations. Alter- 
natively, the large linewidths might indicate that the cores 
are collapsing. In Sect. 14.21 we discuss in more detail the dy- 
namics of the cores in our sample, and relate it to previous 
single-dish studies. 



4 DISCUSSION 

4.1 Temperature of the gas and IR radiation 

In Sect. 13.31 we reported a correlation between the tempera- 
ture of the cores and the luminosity of the region (see Fig. [3] 
panel A). These higher temperatures can result from vari- 
ous heating processes: (i) feedback from low-mass stars, e. g., 
molecular outflows, as found in some cores of AFGL 5142 o r 
20293+5932 fe.g.. lZhang et al.1 12002b IPalau et al.l l2Q07bh : 
(ii) shocks from co nvergent flows (e. g. JCsengeri et al .11 20 Hal ; 
iBeltran et al.ll2012r ); (in ) feedback from high-mass stars, e. g., 
strong radiation (e. g., iLi et al]|2003r ). In the following, we 
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Figure 3. Scatter plots showing the correlations with respect to the bolometric luminosity of (a) the rotational temperature, (b) the 
NH3 (1,1) linewidth, and (c) the core mass, in the top panels; and the correlations with respect to the linewidth of (d) the rotational 
temperature, (e) the ammonia column density, and (f) the core diametre. 'Quiescent starless' and 'perturbed starless' cores are shown 
in red filled and open circles respectively; while 'protostellar' and 'tentatively protostellar' cores are shown as blue filled and open stars, 
respectively (see Sect. 13.11 for details). The black dotted lines show the result of least-squares fits to the data (see Sect. 13. 3|) . with the 
Pearson correlation coefficient (p) shown in the top right corner, together with Kendall's r correlation coefficient. 
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Figure 4. Received flux for each ammonia core versus the ro- 
tational temperature (see Sect. 14.11 for details). The dotted line 
shows a least-squares fit to the data, and the solid line corre- 
sponds to the relation L bol /(47r(A M iR-NH 3 ) 2 ) = csb ^ r 4 ot , with 
ergs the Stephan-Boltzmann constant. Symbols are as in Fig. [3] 



test the feedback from massive stars as a possible agent for 
the increase of temperature in the dense cores. 

Our interferometric observations permit the study of 
the spatial distribution of the cores, and their physical prop- 
erties, with respect to, e. g., the (strong) infrared sources and 
hence the massive stars. In our regions, the bolometric lu- 
minosity is generally dominated, at mid-infrared (MIR) and 
far-infrared (FIR) wavelengths, by a single strong source as- 
sociated with the most massive star. We assume that the 
main source of radiation that can heat the surrounding gas 
is the massive star detectable in the MIR. To determine if 
there is a relation between the temperature of the cores and 
the proximity of the MIR source, w e searched the MSX Point 
Source Catalog ([Price et al.lll999[ ) for the coordinates of the 
MIR sources and calculated the MSX to NH3 core separa- 
tion (Amir-nh 3 ; using the coordinates listed in Table [3] for 
our cores). Using the luminosity and position of the most 
massive (radiating dominant) source and the separation to 
each ammonia dense core, we derived the incident flux re- 
ceived by each core as Lboi/(47r(AMiR-NH 3 ) 2 )- Note that we 
are measuring a projected separation, which is a lower limit 
to the true separation. Thus, for all cores, we actually have 
an upper limit on the received flux, which affects the abso- 
lute normalization of the derived relation and increases the 
scatter of the data points in this relation. In Fig. 2] we plot 
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Figure 5. Scatter plot showing the relatio n between the viria l 
mass M vir (derived as M vir = 210 R Av 2 : iMacLaren et al. 1988) 
and the mass of the core M core . Blue-star and red-circle symbols 
correspond to the protostellar and starless cores studied in this 
work (symbols as in Fig. [3j. Grey squares correspond to single- 
dish ammonia data (Sepulveda et al. 2011) and grey crosses cor- 
respond to 13 CO single-dish data ( Roman-Duval et al.l l2010). The 
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Figur e 6. Var iation of the linewidth, Av, with core radius, fol- 
lowing [CsengerreT^l] (|2011bl) . Red and blue symbols correspond 
to the VLA ammonia data shown in this paper (see Fig. [3] for 
details), grey squares c orrespond to single-dish ammonia data 
(Sepulve da et allfeoill). and grey crosses c orrespond to 13 CO 
single-dish data (JRoman-Duval et al.l l2010h . The dashed line 
indicates a Larson like linewidth-size relation (Av oc size 0-38 ). 
The thick soli d line correspond to the turbulent regulated quasi- 
static model (McKee & Tan 2003), with a surface mass density of 
1.3 g cm -2 and a core mass of 0.1 to 100 Mq (each interval cor- 
responds to 0.1, 1, 10, and 100 Mq). The black cross marks the 
mean and dispersion values of the linewidth (1.36 zb 0.93 km s _1 ) 
and radius (0.032 ± 0.019 pc) for the cores studied in this work. 



the received flux versus the rotational temperature. If the 
radiation of the most massive star dominates the heating 
of its surroundings, we expect that the cores closer to the 
most massive source, i. e., with a large received flux, will have 
higher temperatures. This is consistent with the trend shown 
in Fig. [H where the highest temperatures correspond to the 
cores with a higher received flux. The black solid line cor- 
responds to the relation Lboi/(47r(A M iR-NH 3 ) 2 ) = o"SBT r 4 ot , 
with <7sb the Stephan-Boltzmann constant, while the dotted 
line corresponds to a least-squares fit with a correlation co- 
efficient of 0.6. For a given bolometric luminosity, the tem- 
perature is a factor of ~ 3 times lower for a core located 
at a distance ~ 0.5 pc, than for a similar core located at 
< 0.05 pc. 

Summarising, the radiation deposited in the environ- 
ment by the most luminous star of the cluster can be re- 
sponsible for the temperatures measured in the cores of our 
clustered regions. Thus, the large temperatures measured for 
our starless cores, in comparison with the temperatures mea- 
sured in low-mass star-forming regions, can be interpreted 
as external heating by the radiation of the most massive 
(forming) stars of the cluster. 



4.2 Dynamical state of the cores 

We estimated the virial mass of each core as M v i r = 
210 [ fl/pcl [At; /km s" 1 ] 2 M ^ (|MacLaren et al.1 Il988l : see 
also lBertoldi &; McKed 1992), where the numerical value ac- 
counts for an uniform density distribution across the core, 
and R and Av are the radius of the core and the NH3 (1,1) 
linewidth, respectively (both listed in Table [3]). Note that 
this estimate considers the simplest case of virial equilib- 
rium where only gravity and velocity dispersion of the gas 



are taken into account, thus neglecting the effects of, e.g., 
external pressure and magnetic fields. The calculated virial 
masses are listed in Col. 10 of Tabled In Figure[5] we present 
a plot of virial mass, M v i r , against the mass of the core de- 
rived from the ammonia emission, M core . In Fig. [5] we have 
also included the ammonia clumps studied bv Sepulveda et 
al. (2011) with the Haystack 37 m single-dish telescop e (grey 
squares), and the 13 CO clouds of iRoman-Duval et al.l (|2010h 
observed with the FCRAO 14 m single-dish telescope (grey 
crosses) . These two samples permit to trace clumps at scales 
~ 0.1-10 pc. From the plot we find that the trend outlined 
by this population of clumps is nicely followed by our NH3 
core population sampling lower masses, with the slope of 
the trend close to 1. This, in general, favors a scenario of 
a self-similar, self- gravitating, virialised hierarchy of struc- 
tures from clumps to cores. The fact that the virial masses of 
almost all the cores at 0.05 pc (studied in this work) are sev- 
eral times larger than the core masses (M v i r /M C ore ~ 10Jj 



6 Note that the offset seen in Fig. [5]can not be easily explained in 
terms of the uncertainties of the observational quantities involved 
in the plot. The two main sources of uncertainty are the ammonia 
abundance and the linewidth. For the ammonia abundance we 
assumed a mean value of (4.2 ±2. 7) x 10~ 8 - This value is o btained 
from different survevs (Pillai et al. 2006; Foster et al. 2009; Friesen 
et al. 2009; Rvgl et al. 2010; Chir a et al.ll2013h . over a total of 
122 clumps of different star forming regions. Thus, we expect an 
uncertainty of ~ 60% for the mass of the core. The linewidth is 
obtained from the spectral line fit described in Sect. 13.21 which 
provides errors of a few tenths of km s _1 , corresponding to a 10— 
30% uncertainty in linewidth. These observational errors alone 
are expected to affect the data randomly instead of introducing a 
systematic offset, and thus are not likely the cause of the observed 
offset of about one order of magnitude. 
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suggests that the initial assumptions in our virial analy- 
sis might be too simplistic and that other mechanisms of 
support should be taken into account (e.g., external pres- 
sure, magnetic fields, large scale velocity gradients; see Mac 
Low & Klessen 2004; McKe e fe Ostrikerll 20071 and references 
therein). The turbulence-reg ulated quasi-static scenario pro- 
posed bv lMcKee fe Tanl (J200l l2003h considers the effects of 
the external pressure and an internal pressure that is de- 
rived from a mixture of turbulence and large-scale magnetic 
fields, to characterise the evolution of dense cores until the 
formation of massive stars. Thus, a virialised core that is 
in pressure equilibrium with its natal self- gravitating clump 
has a mean velocity dispersion, cr v ir, and a core radius, i? v ir, 
defined by 



<T v ir 


= 1.09 


iVl core 

60 M 


1/4 


^clump 

g cm -2 


km s _1 _ 



Ry 



pc 



0.057 



M C( 



60 M G 







1/2 



-'clump 



g cm- 



■11/4 



-1/2 



(1) 



(2) 



where the num e rical values appl y to the fiducial cores of 
iMcKee fc Tanl (J2003h (see also iTan et alJ l2013h with a 
power- law density distribution p oc r~ kp with k p = 1.5. 
M core corresponds to the mass of the core, and E c ium P is the 
mass surface density of the clump where the core is embed- 
ded, assumed to be equal to the external pressure ( McKee 
& Tan 2003). In Fig. l6l we follow the a pproach of Csen- 
geri et al. ((2011b|) to compare our observati ons (l inewidth 
and co re size) to the turbulent core model of Mc Kee &; Tanl 
(|2003r ) (represented by the thick black line). The turbulent 
dispersion and the size of the cores of the model have been 
computed for a 0.1 to 100 Mq core mass, and adopting 
Sciump = 1.2 g cm -2 as an average mass surface density of 
the clump, which has been estimated for our regions from 
single-dish (sub) millimetre continuum observations (see Ta- 
ble [6] for details). From Fig.[6l we find that the cores studied 
at scales of 0.05 pc have a mean and dispersion linewidth 
and radius of 1.36=b0.93 km s _1 and 0.032±0.019 pc, respec- 
tively. These values, indicated by a black cross in Fig. [6l are 
slightly below the predictions of the turbulent core m odel, 
a result recently found also bv lCsengeri et al.l (|2011bh in a 
smaller sample of cores in the Cygnus-X region. This sug- 
gests that, when taking into account the external pressure 
effect at scales ~ 0.05 pc, either other forms of support (e. g., 
strong magnetic fields at core scales) are required to stabilise 
the cores, or the cores are close to collapse or even collapsing. 
Observations to gather information on the magnetic field as 
well as on the internal kinematics of the cores should defini- 
tively confirm this picture of the NH3 dense cores at 0.05 pc 
in our sample. 



5 CONCLUSIONS 

We present the results of VLA observations of the NH3 (1,1) 
and (2,2) inversion transition lines towards a sample of 15 
massive star- forming regions. We obtained an average spa- 
tial resolution of ~ 0.05 pc, and analysed the ammonia emis- 
sion in a consistent way. Our main conclusions can be sum- 
marised as follows: 



Table 6. Parameters derived for the clumps associated with the 
15 massive star-forming regions 



S u h Mass c FWHM d 



# a (Jy) (M ) 



n 



diametre d 

(") / (p<0 



^clump 

(g cm" 2 ) 



1. 


1.2 


85 


32.4 


30.5 


/ 


0.27 


0.32 


2. 


15.3 


278 


28.7 


25.1 


/ 


0.22 


1.55 


3. 


6.8 


123 


44.2 


41.9 


/ 


0.37 


0.25 


4. 


21.1 


384 


40.2 


37.7 


/ 


0.33 


0.94 


5. 


4.5 


64 


30.1 


26.6 


/ 


0.21 


0.40 


6. 


7.3 


199 


31.9 


28.7 


/ 


0.31 


0.57 


7. 


32.9 


1160 


30.0 


26.5 


/ 


0.32 


2.97 


8. 


13.3 


201 


34.5 


31.5 


/ 


0.25 


0.85 


9. 


28.3 


2290 


37.8 


35.1 


/ 


0.65 


1.46 


10. 


15.2 


341 


45.0 


42.8 


/ 


0.42 


0.53 


11. 


15.8 


174 


56.5 


54.7 


/ 


0.37 


0.34 


12. 


6.7 


253 


38.5 


35.9 


/ 


0.45 


0.33 


13. 


4.5 


145 


34.6 


31.6 


/ 


0.37 


0.29 


14. 


5.2 


17 


16.5 


8.8 


/ 


0.03 


4.30 


15. 


67.1 


185 


41.6 


39.2 


/ 


0.13 


2.78 



a Number of the region as indicated in Table [T] 

b Flux d ensity at 850 ^m ob tained from the SCUBA Legacy 

Project (|Di Francesco et aljfiooi ), exc ept for: 00117+6412 

(region #1: 1.2 mm from IRAM 30 m, Sanchez-Monge et al. 

2008 ): 20081+3122 (region #7: 850 urn from JCMT. Thompson 

et ahl2006h; 22198+ 6336 (region #14: 850 /im from JCMT, 

Jenn ess et al.l ll995) . 

c Mass of the clump measured from the JCMT at 850 fim or the 

IRAM 30 m at 1.2 mm, assuming a dust temperature of 20 K, a 

dust (+gas) mass opacity coefficient of .0175 cm 2 g -1 at 



Hennina 



850 | im and 0.0090 cm 2 g _1 at 1.2 mm (jOssenkopf . 
Il994l) , and a gas-to-dust mass ratio of 100. 
d Observed diametre at full width half maximum obtained 
(FWHM) from a Gaussian fit, and deconvolved diametre derived 
taking into account a main beam for the SCUBA observations 
at 850 /im of 14". Our sizes obtained from a Gaussian f it are 
slightly different to those reported bv lPalau et al.l (|2013T) , which 
were estimated by measuring the observed diametre directly on 
the image, and using a beam size of 23" for deconvolution, 
which affects mainly the sources with the smallest sizes. 
e Clump mass surface density derived as S c lump = Mass/Area 
with Area = ir (diametre/2) 2 . 



- We identified a total of 73 ammonia dense cores us- 
ing the two-dimensional CLUMPFIND algorithm. On av- 
erage we found of the order of 5-6 cores in each field (of 
~ 2 arcmin, or 1.5 pc at a distance of 2.5 kpc), with sizes 
~ 0.06 pc. The cores were classified as quiescent starless, 
perturbed starless, and protostellar based on its association 
with infrared emission and taking into account other star 
formation signposts (e.g., outflows, masers, centimetre con- 
tinuum emission) reported in the literature. 

- The average values of linewidth, opacity, rotational 
temperature, and ammonia column density for the 73 dense 
cores are ~ 1.4 km s _1 , 1.4, 19 K, and 10 15 cm -2 , respec- 
tively. We found that protostellar cores have larger temper- 
atures and linewidths (~ 21.3 K and 1.8 km s _1 ) than qui- 
escent starless cores (~ 16.0 K and 1.0 km s _1 ). Similarly, 
cores classified as perturbed starless have temperatures and 
linewidths (~ 19.3 K and 1.4 km s _1 ) slightly larger than 
quiescent starless cores, but smaller than cores harboring a 
protostar. 
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- The average temperatures derived for the quiescent 
starless cores in our clustered regions (~ 16 K) are larger 
than those measured in starless cores in more isolated re- 
gions (~ 10 K). Our starless and protostellar cores show 
important non-thermal supersonic components (cr n th ~ 
0.44 km s _1 and 0.76 km s _1 , respectively) and cr n th/c s ~ 
2.0, well above to the close-to-thermal lines measured in iso- 
lated star- forming regions. 

- We found that the temperature and NH3 (1,1) linewidth 
of the cores are correlated with the luminosity of the region, 
suggesting that the feedback of the luminous massive (form- 
ing) stars can affect the surrounding molecular environment. 

- We evaluated a possible relation between the tempera- 
ture of the core and its separation with the most important 
radiating source in the region, in terms of the incident flux 
received by the core. The highest temperatures are mea- 
sured for those cores located close to the brightest MIR 
source (separations < 0.05 pc), with the temperature de- 
caying by a factor of ~ 3 when the cores are located at 
distances ~ 0.5 pc. Thus, the large temperatures found in 
starless cores in massive star-forming regions, can be in- 
terpreted as external heating by the radiation of most the 
massive star of the cluster. 

- We compared the virial masses of our cores (at scales 
~ 0.05 pc) with those of clumps with sizes ~ 0.1-10 pc. A 
simple virial equilibrium analysis seems to suggest a trend 
between M v i r and core/clump masses favoring a scenario of 
a self-similar, self-graviting, virialised hierarchy of structures 
from clumps to cores. A closer inspection of the dynamical 
state at core scales, taking into account external pressure 
([McKee &; Tad 120031 ), indicates that either another form of 
support (e.g., from magnetic field) is required to stabilise 
the cores, or that the cores are close to collapse or even 
collapsing. 

Overall, our data suggest that the cores forming in clus- 
tered environments harboring intermediate/high-mass stars 
are more affected/perturbed than cores in clustered low- 
mass star-forming regions, and thus the initial conditions for 
star formation in these environments seem to be different to 
those typically assumed in low-mass star- forming regions. 
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APPENDIX A: HYPERFINE STRUCTURE 
(HFS) FITTING 

Al Parameters of the fit 

The usual fitting of the hyperfine structure of a transition, 
for instance, the inversion transition NH3 (J,K) — (1,1), 
estimates four free parameters of the transition, namely 

i) r m , the optical depth of the main component, 
ii) A x r m , product of the amplitude (A = f[J v (T e ^) — 
Jzy(Tbg)], where / is the filling factor, J»(T) is the Planck 
function in units of temperature, T ex is the excitation tem- 
perature, and Tb g is the background temperature) times the 
optical depth of the main component, 

iii) ^lsr, the central velocity of the main component, and 
iv) Av, the linewidth of the hyperfine components, 

with the usual assumptions that T ex , vlsr, and Av are the 
same for all the hyperfine components. 

However, some of the parameters are ill-conditioned in 
some cases. In the optically thin limit, r m < 1, only an up- 
per limit for the value of r m can be obtained. However, the 
parameter A x T m gives the peak intensity of the line (pro- 
vided that the hyperfine components are unresolved, or, in 
the case NH3 , the hyperfine magnetic components that form 
part of the main component), and is well constrained. This 
parameter, in the Ray leigh- Jeans approximation, is neces- 
sary to calculate the column density of emitting molecules. 
In the optically thick limit, r m ^> 1, only lower limits for 
r m and A x r m can be obtained (and consequently, a lower 
limit for the column density). 

A better fit parameter than T m , with well defined values 
in the optically thin and optically thick limits, is 1 — e~ Trn . 
This parameter always has values between (optically thin 
limit, Tm <C 1) and 1 (optically thick limit, T m ^> 1)- Sim- 
ilarly, A (l — e _Trn ) is a better fit parameter than A x r m , 
since it is well constrained in all cases. 

The Hyperfine Struct urqj (hereafter HfS) fitting proce- 
dure used to fit the NH3 (J, K) — (1, 1) lines in this paper, 
fits the four parameters 

i) 1 - e~ r ™, 
ii) A(l-e- T "), 
iii) ^lsr, 
iv) Av. 



Table Al. Values of A(ra, a) for calculating the parameter un- 
certainties, where m is the number of fitted parameters simulta- 
neously, and a is the significance level, given in percent and in the 
equivalent number of sigmas for a Gaussian error distribution. 







a 








68.23% 


90.00% 


99.00% 


99.90% 


m 


(la) 


(1.64<t) 


(2.57cr) 


(3.30cr) 


1 


1.00 


2.71 


6.64 


10.83 


2 


2.30 


4.60 


9.21 


13.82 


3 


3.53 


6.25 


11.34 


16.27 


4 


4.72 


7.78 


13.28 


18.46 


5 


5.89 


9.24 


15.09 


20.52 


6 


7.04 


10.64 


16.81 


22.46 


7 


8.17 


12.02 


18.48 


24.32 


8 


9.30 


13.36 


20.09 


26.12 


9 


10.42 


14.68 


21.67 


27.87 


10 


11.53 


15.99 


23.21 


29.59 



The values and uncertainties for the parameters needed to 
derive physical properties of the emitting gas, r m and A x 
T m (or their limits in the two limiting cases) can be easily 
derived from the values and uncertainties for the parameters 
of the fit. 



A2 Fitting procedure 

T he fitting strategy was similar to that used by Estalella et 
al. (|2012i ). HfS samples the parameter space of dimension 
four, defined by the parameters 1 — e~ Trn , A (l — e _T ™), 
^lsr, and Av, to find the minimum value of x 2 , 



E 



y?°\pi 



(Al) 



where yf° s are the observed line intensities for the n spectral 
channels, yT° d {Pi, • • • )Pm) are the model line intensities, de- 
pending on m free parameters pk, and <ii are the errors of 
the observations. 

If the model is a good fit to the observations, x 2 has a 
minimum value given by 



2 

Amin 



(A2) 



This can be used to estimate the observational error, a (the 
same for all points), if unknown, 



1 n 

— Ek 

— m ^-^ 1 



Vi 



(A3) 



Several sampling methods of the m-dimensional pa- 
rameter space are possible, i. e., regular grid, random, Hal- 
ton sequence (JHaltonJ 1964). We adopted a Halt on quasi- 
random sequence because it samples the parameter space 
more evenly than a purely random sequence, and the con- 
vergence of the fitting procedure to the minimum of \ 2 is 
faster. 



A3 Parameter error estimation 



7 HfS is freely available and can be downloaded 
http://www.am.ub.edu/~robert/hfs/hfs.html 



from 



Once the best fit parameter values for which % 2 

are found, the uncertainty in the fitted parameters, a(pk 



Amin 
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(k = 1 . . . ?n ), can be estimated as the projection over each 
parameter axis of the confidence region of the parameters, 
the region of the m- dimensional parameter space for which 
X 2 does not exceed the minimum value by a certain amount, 

Xmin + A. 

Following IXvn3 (|l976h and lWall fc Jenkinsl (J2003I ). the 
probability 

Prob [x 2 - Xmi„ sS A(m, a)] = a, (A4) 

is that of a chi-square distribution with m degrees of free- 
dom, and a is the significance level (0 < a < 1). Thus, 
A(ra, a) is the increment in % 2 such that if the observation 
is repeated a large number of times, a fraction a of times 
the values of the fitted parameters will be inside the confi- 
dence region, i.e. in the interval pk =b cr(pk)- For instance, in 
the case of this paper, m — 4, and for a significance level of 
0.68 (equivalent to 1 a for a Gaussian error distribution), the 
valu e of the increment in viL„ is A (4, 0.68 ) = 4.72 (X ampton 
et al. ll976l : see Table [XT]) 

Assuming that the model is a good fit to the observa- 
tions, i. e., Xmin — n ~ m i the condition 

X 2 = Xmin + A(m,oO, (A5) 

-f ,1+^^. (A6) 

X m in ri-m 

This expression is useful since it can be given in terms of 
the rms fit residual, e, 

2 _ E: = i[(^ bs -^ mod )M] 2 _ x 2 , A ~ 

eiuim 2 _ e: =1 im 2 ' ( } 

for which we obtain that the confidence region is given by 
the parameter values that increase the rms fit residual to 



can be written as 



1+ AKa) 

n — m 

This last expression can be used even when the errors of the 
observations are unknown. 



APPENDIX B: VELOCITY, LINEWIDTH AND 
SPECTRA IMAGES 



In Table IBlt we summarize the infrared properties of each 
core, together with the information of other star formation 
signposts reported in the literature. In Col. 4 of the table, 
we classify each core as 'quiescent starless' (S), 'perturbed 
starless' (S*), 'protostellar' ( P) or 'tentatively protostellar' 

(p?). 

In Fig. lBU we present the first-order (intensity weighted 
mean vlsr) and second-order (intensity weighted mean 
linewidth) moment maps for the 15 massive star- forming 
regions studied in this work. Note that the second-order 
moment maps have been corrected to show directly the 
linewidth (Av) instead of the dispersion (a v ), by using the 
relation Av — VS In 2 a v . 

In Fig. IB21 we present the observed and fitted spectra 
of the NH3 (1,1) and NH3 (2,2) transitions for all the cores 
identified and listed in Table [3] 
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Table Bl. Classification of ammonia cores in protostellar or starless 



# Type b IRAC MIPS/WISE mm cm outflow maser comments 



ID Region a 



01 00117+6412 


01 


P 


point 


- 


yes 


no 


no 


no 


02 


02 


S 


dark 


- 


yes 


no 


no 


no 


03 


03 


P 


point 


yes 


yes 


yes 


yes 


yes 


04 


04 


s* 


dark 


- 


yes 


no 


no 


yes 


05 


05 


p? 


point 


- 


no 


no 


no 


no 


06 AFGL5142 


01 


p 


point 


yes 


yes 


yes 


yes 


yes 


07 


02 


p? 


point 


- 


no 


no 


OA 


no 


08 


03 


p? 


point 


- 


no 


no 


no 


no 


09 05345+3157NE 


01 


p 


point 


- 


yes 


yes 


no? 


no 


10 


02 


p 


point 


yes 


yes 


no 


yes 


yes 


11 


03 


s 


dark 


- 


no 


no 


no 


no 


12 


04 


s 


dark? 


yes? 


yes 


no 


no 


no 


13 


05 


s 


dark 


- 


yes 


no 


no 


no 


14 


06 


s 


dark 


- 


yes 


no 


no 


no 


15 


07 


s 


dark 


- 


no 


no 


no 


no 


16 05358+3543NE 


01 


p 


ext? 


- 


yes 


no 


yes 


no 


17 


02 


p? 


point 


- 


no 


no 


OA 


no 


18 


03 


p 


point 


yes 


yes 


yes 


yes 


yes 


19 


04 


s* 


dark 


- 


no 


no 


OA 


no 


20 


05 


p 


point 


- 


no 


no 


yes 


no 


21 


06 


s* 


dark 


- 


no 


no 


OA 


no 


22 


07 


s 


dark 


- 


no 


no 


no 


no 


23 05373+2349 


01 


p 


point 


- 


yes 


yes 


yes 


no 


25 


02 


s* 


dark? 


yes? 


faint 


no 


no 


no 


27 19035+0641 


01 


p 


point 


yes 




yes 






28 


02 


s* 


dark 


- 




no 






29 


03 


s 


dark 


- 




no 






30 


04 


s? 


ext 


- 




no 






31 


05 


s 


dark 


- 




no 






32 20081+3122 


01 


p? 


point 


- 


yes? 


no 


OA 


no 


33 


02 


s 


dark 


- 


no 


no 


no 


no 


35 20126+4104 


01 


p 


point 


yes 


yes 


yes 


yes 


yes 


36 


02 


s 


dark 


- 


yes 


no 


no 


no 


37 


03 


p? 


point 


- 


no 


no 


no 


no 


38 G75.78+0.34 


01 


p 


point 


yes 


yes 


yes 


yes 


yes 


39 


02 


s* 


ext 


- 


no 


no 


no 


no 


40 


03 


p? 


point 


- 


no 


no 


no 


no 


41 


04 


s 


dark 


- 


no 


no 


no 


no 


43 


05 


s 


dark 


- 


no 


no 


no 


no 


44 


06 


p? 


point 


- 


no 


no 


no 


no 


45 20293+3952 


01 


s* 


dark? 


- 


no 


no 


no 


no 


46 


02 


s* 


dark? 


- 


yes 


no 


OA 


no 


47 


03 


p 


point 


yes 


yes 


no 


yes 


yes 


48 


04 


p? 


point 


- 


yes 


no 


no 


no 


49 


05 


p 


point 


yes 


yes 


no 


yes 


no 


50 


06 


p? 


point 


- 




no 


no 




51 20343+4129 


01 


s* 


ext 


- 


yes 


no 


no 


no 


52 


02 


s* 


dark 


- 


yes 


no 


OA? 


no 


53 


03 


p? 


dark? 


- 


yes 


yes? 


OA? 


no 


54 


04 


p 


point 


yes 


faint 


yes 


no 


no 


55 22134+5834 


01 


s 


dark 


- 




no 




no 


56 


02 


s* 


dark? 


- 




no 




no 


57 


03 


s 


dark 


- 




no 




no 


58 


04 


s 


dark 


- 




no 




no 


59 


05 


p? 


point 


- 




no 




no 



no mm/outflow/maser associated 
core along outflow axis 

point IR source at the border, but dark 

core along outflow axis 
near CO outflow lobe 
near SiO outflow lobe 

between two strong IR sources 
close to bright IRAC source 

mm source A" to the south of the polygon 



near strong IRAC source 



near strong IRAC source 
perturbed by outflow 



near strong IRAC source, perturbed? 
near strong IRAC source, perturbed? 
near strong IRAC source, perturbed? 



near strong IRAC source, perturbed? 



20 A. Sanchez- Monge et al. 



Table Bl - continued 



ID Region a 


# 


Type b 


IRAC 


MIPS/WISE 


mm 


cm 


outflow 


maser 


comments 


60 22172+5549N 


01 


P 


point 


yes 


yes 


no 


yes 


no 




61 


02 


S 


dark 


- 


no 


no 


no 


no 




62 


03 


P? 


point 


yes 


faint 


no 


OA 


no 


core along outflow axis 


63 22198+6336 


01 


s 


dark 


- 


no 


no 


no 


no 




64 


02 


p? 


point 


- 


no 


no 


no 


no 




65 


03 


p 


point 


yes 


yes 


yes 


yes 


yes 




66 


04 


s* 


ext 


- 


no 


no 


OA 


no 


core along outflow axis 


67 


05 


s 


dark 


- 


no 


no 


no 


no 




68 


06 


s 


dark 


- 


no 


no 


no 


no 




69 


07 


s 


dark 


- 


yes 


no 


no 


no 




70 


08 


s 


dark 


- 




no 




no 




71 CepA 


01 


s* 


dark 


- 




no 


OA 


no 


core close to outflow axis 


72 


02 


s* 


dark 


- 




no 


no 


no 




73 


03 


s* 


dark 


- 




no 


OA 


no 


core close to outflow axis 


74 


04 


p? 


point 


- 




no 


OA 


no 


core close to outflow axis 


75 


05 


s* 


ext 


- 




no 


OA 


no 
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a References: We have made use of the Spitzer /IRAC, Spitzer /MIPS, and WISE images to characterise the IR properties of the cores 

at wavelengths ranging from 3.6 fim to 24 jum. For each region we have searched the literature for information on millimeter 

continuum, centimetre continuum, outflow and maser emission. The re ferences used (for each region) are: 

00117+6412 : IPalau et all (2010); Busquetl (l2010h ; TSanchez-Mong j (l201ll); 

AFGL5142: Zha ng et al. | J2 002. 2007); Busaue t](|2010h : ISanchez-Mongd d201lh : iBusauet et alJ (EoTlh : IPalau et alJ (EoTll . l2013h : 

05345+3157NE: Fon tani et aLN 2008. 2009, 2012a)j 

05358+3543NE: Beut her et alJ d2002cL l2007h: iLeurini et al.1 (l2007h: ISanchez-Mongel d201lh: 

0537 3+2349: iRodriguez et al.l dl980h: IVarricatt et al.l (l2oToh: iMolinari et al.l dl996U2002U2008h: iGutermuth et alJ d2009h: Khanzadvan 

et al. (|201lh : 

19035+0641: 

20081+3122: 

20126+4104: 



Sanchez-Mongd J201ll) 



Su et al 



al. 
(2. 



( 2009h:lKumar et al.l l2004 i; 
< '• .w. ni et al.ldl997l. 11 999. 2005) ~lKeto & Zhang! fcOlCh : IBusauetl d2010h : 



G75.78+0.34 ISanchez-Mongd d201lh: ISanchez-Monge et all d2013h : 
20293+3952: iBeuther et al. (2002d, 2004ab); Palau et al.l d2007al) 



20343+4129: Palau et al.1 I <2Q07bh: fr ont ani e t alJ J2qi2h 

22134+5834: IBusauetl J2010h : I Sanchez-Mongd d201lh: 

22172+5549N: Font ani et atl (120041): IPalau et al.l d201ll. 12013^: 

22198+6336: Sanchez-Monge et al. (2010): ISanchez-Mongd (|201lh : IPalau et al.1 (J201lil20ld) : 

CepA: iTorrelles et al.l dl986h : iPatel et al.1 (|2005h . 

b Type of the core: S: quiescent starless, S*: perturbed starless, P: protostellar, P?: tentatively protostellar (see Sect. 13. l|) . 
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Figure Bl. Velocity and linewidth maps (in units of km s — 1 ) obtained from the first and second-order moment maps. Note that the 
second-order moment map has been multiplied by V8 In 2 to obtain linewidth instead of dispersion. For each region we show two panels 
corresponding to the intensity weighted mean ^lsr (top panel) and to the intensity weighted mean velocity linewidth (bottom panel). 
Colour scales are indicated in the top of each panel. Symbols are as in Fig. [T] 
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Figure B2. Observed (black histogram) and fitted (thick blue line) spectra of the NH3 (1,1) and (2,2) lines for the cores listed in TableE] 
Vertical axis corresponds to brightness temperature (in K), and horizontal axis corresponds to velocity (in km s — 1 ). 
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